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ABSTRACT 

We review the progress in understanding the strange particle yields in nuclear colli- 
sions and their role in signalling quark-gluon plasma formation. We report on new 
insights into the formation mechanisms of strange particles during ultrarelativistic 
heavy-ion collisions and discuss interesting new details of the strangeness phase di- 
agram. In the main part of the review we show how the measured (multi-)strange 
particle abundances can be used as a testing ground for chemical equilibration in 
nuclear collisions, and how the results of such an analysis lead to important con- 
straints on the collision dynamics and space-time evolution of high energy heavy-ion 
reactions. 

1 Introduction 

The high energy heavy-ion program at the Brookhaven Alternating Gradient Syn- 
chrotron (AGS) and at the CERN Super Proton Synchrotron (SPS) is devoted to the 
investigation of nuclear matter under extreme conditions (high temperatures and/or 
densities). The most exciting prospect is the possible observation of a phase tran- 
sition from normal nuclear matter to a so-called Quark-Gluon Plasma^ (QGP). A 
large amount of experimental and theoretical effort has gone into working out clear 
signatures for the QGP. Although the ultimate hope is to eventually see the QGP 
"shine" via direct electromagnetic radiation of photons and lepton pairs^, for prac- 
tical reasons (larger signals and better understanding of backgrounds) experimental 
efforts have up to now focussed on the production of hadrons with newly produced 
quark flavors (strangeness, charm, etc.). At present the largest body of data with the 
so far most promising indications for unconventional physics has been accumulated 
in the sector of strange particles, following the early suggestions of Rafelski^ well over 
10 years ago. 

There are two major aspects of the so-called "strangeness signal": The more spec- 
ulative one is connected with the search for exotic strange matter, which requires the 
intermediate production of a QGP state. The list of these so far undetected phenom- 
ena contains strangelets'*"^ (small droplets of metastable strange quark matter - the 
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ultimate proof of the existence the quark-gluon plasma), but also the H-dibaryon 
and more recently^ the so-called MEMOs (Metastable Exotic Multistrange Objects). 
For this aspect of the strangeness signal we refer the reader to Refs. 9, 10 and to the 
contribution by C. Greiner and J. Schaffner in this volume. 

The other, more conventional aspect of strangeness as a QGP signature was in- 
troduced by Rafelski who predicted an enhanced production of strange hadrons in 
heavy-ion collisions when compared to nucleon-nucleon or nucleon-nucleus collisions 
at the same energy^'^^. The idea is simple and is based on the different produc- 
tion mechanisms of strangeness in a QGP (individual strange quark pair production 
from a dense system of gluons and light quarks) and in a hadron gas (production 
of hadrons with opposite strangeness in inelastic hadron-hadron interactions), which 
have strongly different intrinsic time scales. But the quantitative interpretation of the 
observed enhanced nuclear strangeness production as a QGP signal is a subtle issue. 
Clear and unambiguous statements require a quantitative theoretical understanding 
of strangeness production in a QGP and in a hadronic environment (including the 
initial pre-equilibrium stages), of the evolution of strangeness during the fireball life 
time, in particular during the difficult hadronization stage, and of the variations of 
the non-strange multiplicities as one goes from nucleon-nucleon to nuclear collisions. 
Some of these issues are still not satisfactorily resolved, but significant progress has 
been made in the last few years which we will review in this article. 

An important, not always appreciated advantage of strangeness is the large variety 
of strange particles species. Due to the long life time of strange quarks against weak 
decay, one has by now been able to identify and measure the yields and spectra of 
10 different species of strange hadrons in heavy-ion experiments (if"*", K~ , K°, ^, A, 
A, S~, S"'", Q'^). An up-to-date list of references is given in Tables 8 and 9 in 
the Appendix; an account of the latest results can be found in Refs. 12, 13. These 
data triggered a lot of studies within various models for nuclear collisions. Especially 
thermal models^'* which are characterized by a rather small number of adjustable 
parameters and thus predict various restrictive relations between the various particle 
species can thus be tested very precisely. We will devote a large fraction of this article 
to these investigations, not only because large progress has been made in particular 
in this sector, but also because these studies provide a simple and intuitive picture 
which gives valuable qualitative and quantitative insight into the physics of hadron 
production in nuclear collisions. The thermal picture allows for a straightforward 
investigation of the question of chemical equilibration in the strange and non-strange 
sectors and provides a simple characterization of the final stages of the the heavy-ion 
collision in terms of freeze-out parameters, entropy production, and collective fiow 
dynamics. It provides an intuitive bridge between the data on the one hand and 
complete kinetic simulations of the phase-space evolution of the collision region on 
the other hand. 

In this review we concentrate on the developments after the publication of the 
first volume of this book-^; for the earlier work we refer to the reviews by Koch^'^^. 
In addition to the discussion of strange hadron yields and momentum spectra we 
also report on some recent progress on the the microscopic mechanisms of strange 
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particle production in a hot and dense environment and on the general structure of 
the strangeness phase diagram. 

The article is divided into the following parts. In Section 2 we cover some new 
developments on strangeness production mechanisms. Section 3 provides the basic 
equations for the thermodynamics of hadronic matter. These are used in Section 4 
to study the phase diagram of strange matter. In Section 5 we discuss additional 
ingredients like collective flow and deviations from chemical equilibrium which are 
necessary when testing the thermal model against experimental data. This is then 
done in Section 6 for various heavy experiments on strange particle production. In 
Section 7 we shortly discuss the aspect of entropy production in these heavy-ion 
experiments. A comprehensive discussion of the outcome of the thermal analysis 
follows in Section 8. In Section 9 we draw some speculative conclusions from this 
discussion on the nature of the hadronization and freeze-out processes and discuss 
the internal consistency of our interpretation of the strange particle abundances. We 
conclude with an outlook in Section 10. 

2 Strangeness Production Mechanisms 

The usefulness of strange particles as probes for the physics of a nuclear collision 
arises from certain characteristic features in their production and kinetic evolution. 
In contrast to pions, for example, which are the most efficient carriers of entropy 
and whose final abundance is thus more or less determined already in the very first, 
hard collision stage of the reaction where most of the entropy is produced^^ - the 
transparency of nuclear matter is the important ingredient -, the abundance of strange 
quarks continues to evolve until the very end of the collision fireball. Thus all stages 
contribute to the final strangeness yield and must be considered. 

The first pre-equilibrium stage is dominated by the physics of hard processes. 
The amount of strange particles created in these hard collisions is rather uncertain. 
For the present highest energy range [^/s = 20 A GeV) this phase is difficult to 
handle quantitatively because it is still dominated by non-perturbative QCD. At 
higher energies the parton model can be more reliably used; at ^/s = 200 A GeV, for 
example, the amount of strangeness produced in the very early collisions has been 
calculated in the Parton Cascade Model to be about 40 % of the equilibrium value^^; in 
the HIJING model, which uses a different soft momentum cutoff procedure, the initial 
strangeness fraction is somewhat lower^^. For the lower AGS and SPS energies one 
must instead make an educated guess, and one generally assumes an initial fraction 
of strange particles as observed in the final state of p+p collisions at the same energy, 
after correcting for resonance decays^^, which corresponds to roughly 20% of the 
equilibrium value^°. 

If the collision region is initially in the deconfined phase, thermal equilibration 
occurs at a very rapid rate^^'^^. Chemical equilibration of gluons and light quarks 
is expected to follow shortly after^-^ although this question is still debated among 
the various high energy event generators^^'^^. Ideally one would thus after about 2-3 
fm/c create a thermally and chemically equilibrated QGP consisting of (mostly) light 
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quarks and gluons. Due to the strange quark mass threshold, strangeness production 
only follows later. Given the above scenario, the production rate and strangeness 
equilibration time scale can be calculated with the methods of thermal equilibrium 
field theory, and if the temperature is sufficiently high, even perturbative QCD can 
be applied. While the first calculations used bare quark and gluon propagators^^ 
recent investigations have begun to incorporate the new knowledge of how to include 
certain collective and non-perturbative effects at high temperature via resummation 
of the perturbative series. These new developments will be described in Section 2.1. 

An alternative scenario is a nuclear collision without a deconfining phase transi- 
tion. Then the relevant degrees of freedom are the hadrons. High temperatures and 
densities will affect the hadron dynamics, their masses and their resonance widths, 
and thus the strangeness production. Of particular relevance is a possible decrease of 
the hadron masses. In Section 2.2 we will report some calculations which take into 
account this effect for strangeness production. 



2.1 Quark-Gluon Production Mechanisms 

The production rates of strange quarks in a QGP are calculated by perturbative QCD. 
At finite temperature, however, any massless theory is harmed by severe infrared 
singularities which are difficult to deal with, especially for gauge theories. Braaten 
and Pisarski^'* developed a resummation technique for the endangered low momentum 
modes at finite temperature and density calculations. The main idea is to use for 
light particles with momenta of order q< gT [g <^ 1) a resummed propagator which 
contains self-energy corrections due to the medium. In lowest order of the g the low 
momentum particles in a thermal environment thereby acquire a thermal mass mth, 
which is given by 

m, = (iVe + Nf/2)g'T'/9 (1) 

for gluons and by 

m, = ml + g'T'/Q (2) 

for quarks. These "massive" quark and gluon modes of finite temperature QCD are 
analogues of the well-known plasmon in the theory of electron gases and electromag- 
netic plasmas. At next to leading order in g these modes also acquire a thermal 
width. In addition to the low-momentum propagators, also some vertices in which 
all momenta are soft of order gT must be resummed^'*. 

The new technique was used to recalculate the production rates for strange quarks 
m a 

QQp35,36 rj.^^ Feynman dia grams for ss production in lowest order QCD are 
shown in Figure 1. There are three different groups: quark anti-quark collisions (a), 
gluon fusion (b) and thermal gluon decay (c). The third process, which is not allowed 
for a bare gluon propagator [nig = 0), is now possible because of the finite mass and 
width of the thermal gluons. 

The expressions for the rates of the processes depicted in Figure 1 are given by^^: 

o _ f d% d% d^ps d^p, . .4c.p IP p p^ 

Ugg^s, - J (27r)32E, {2^Y2E, {2^Y2E, {2^Y2E, ^ ^ ^ " ^ ^ ' 

xMEg)ME,){l - fF{Es)){l - ME,))Y,\M{qq ^ ss)\\ (3) 
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Rgg^s-s = - J (2TrY2E, (2TrY2E, (2TrY2E, ^ ' + P2 - Ps - P.) 



xfB{Er)fB{E2){l - fF{Es)){l - fF{E,)) ^ \M{gg ^ ..-)|^ (4) 



and 



-^g—^ss 



f d^V-s (n^\ian - P - P\ 

J (27r)32E, (27r)32E, (27r)32E, ^ ^ ' 

xfB{E,){l - fF{Es)){l - fF{E,)) E \M{g ^ s-s)\\ (5) 

Capital letters correspond to four- moment a and (/f) denote the Bose- Einstein 
(Fermi-Dirac) distributions, respectively. 

An extensive discussion of the calculation of the first two processes of Figures la) 
and lb) with bare propagators and vertices can be found in the early literature (see 
Refs. 9, 27, 33). A complete new calculation of the relevant matrix elements with 
resummed propagators and vertices is not easy and therefore has so far not been 
presented; work in this direction is in progress^^. The new calculations published 
so far^^~^^ simply replace all bare propagators in Figure 1 by thermal propagators 
with a thermal mass and width. In this approximation the thermal gluon decay. 
Figure Ic), is then the most interesting new phenomenon and has therefore received 
most of the attention. In view of the incompleteness of these calculations it is not 
surprising that the results are controversial. The thermal gluon decay into strange 
quark pairs was first evaluated in Ref. 38. The motivation for a finite gluon mass was 
at that time taken from lattice QCD results^^ rather than from resummed thermal 
perturbation theory. Gluons with masses around 500 MeV fit the data on energy 
density and pressure of a SU(3)c lattice calculation^^ above Tc very well. The rate 
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Figure 2: Partial and total rates 
A = d^Ns/dx'^ for strange quark 
production, according to the calcu- 
lation in Ref. 38: qq —f ss (dashed 
line), gg —f ss (dotted line), g —f ss 
(dashed-dotted line) and M denotes 
the used constant thermal gluon 
mass. 



(5) was calculated with a constant, temperature independent gluon mass, and the 
quarks have been assumed to have their constituent masses. 

The result is shown in Figure 2. The upper plot shows the old calculation^^''*^ with 
nig = (bare gluon propagator), the lower plot a calculation with nig = 500 MeV. In 
the first case the gluon fusion process is by far the dominant source of strange quark 
production, dominating the production from light quarks by an order of magnitude. 
A finite gluon mass decreases the gluon fusion rate, because of the strongly reduced 
Bose-factors (see Eq. (4)) for the incoming gluons. On the other hand it raises the 
strangeness production rate from light qq pairs slightly because the intermediate gluon 
is now less ofF-shell due to its thermal mass. The dominant process, however, is now 
the thermal gluon decay, because (in contrast to the gluon fusion) only one thermal 
distribution function appears in the incoming channel (see Eq. (5)). Note also that 
this process is of order while the other processes are of order g"^. Surprisingly, in 
spite of these drastic changes in the individual production channels, the total rate 

-^tot '^QQ — ^ss -^gg — >ss -^g — >ss 

stays nearly unchanged^^. The authors conclude that the statements about the equi- 
libration time scales in Ref. 33 are nearly unchanged and the gluons remain the main 
strangeness producer. 

Two recent calculations^^'^^ tried to improve this analysis by using the thermal 
propagators of Ref. 34. We will shortly describe the calculation of Ref. 36. The 



thermal gluon propagator may be written as 

iD'liiqo, q) = -iS-' [Pj, Ar(go, q) + Ai(go, q)] , (7) 

where P^ and P^^^ are transverse and longitudinal projectors, respectively, and the 
corresponding propagators are 

Ar,L(go,g) = ^7^-7 r • (8) 

^T,L{qo, q) denotes the self-energy correction due to a thermal loop^'*. The poles in 
the propagator (8) are determined by the dispersion relation 

ql = q' + UTAqo,q)- (9) 

The poles are located at 

go = I^T,L + ilT,L, (10) 

where the imaginary shift of the pole ')t,l is related to the imaginary part of the self 
energy through 

1T,L = Res(Ar,L) Imllr.L , (11) 
with the residue of the propagator at the pole given by 



Res(Ar,L)-' 



dqo 



(12) 



The expression for the gluon propagator can then be rewritten as 



A(Q) = ^^-R^n + ^lRes{Ar .^3. 

(g2_Ren)2 + 72Res(A)-2 ^ (g2_Ren)2 + 72Res(A)-2 ' ^ ^ 

where we suppressed the subscripts T, L. This expression is used to replace the mass- 
shell ^-function for thermal gluons: 

xrD2_^2x 1 7r,LRes(Ar,L)"' . . 

The evaluation of Eq. (5) is then done^^ neglecting the Pauli blocking factors. The 
result is^^'^^ 



Pj_, = T^/ dqofsiqo) dqq'{Q' + 2ml)Jl 

6Tr^ J2m, Jo V 

(15) 



7rRes(Ar)-^ 



(g2 _ Renr)2 +7|Res(Ar)-2 ' 
plus a similar expression for the decay rate of the longitudinal gluons 



7 



The biggest uncertainty in such a calculation is the determination of the imaginary 
part of the pole ')t,l in the gluon propagator. 7^ is the damping rate of transverse 
plasma oscillations in a QGP at high temperature. For plasmons with nonzero spatial 
momentum relative to the heat bath 7^ is not perturbatively calculable, but can be 
related'*-^ to the so-called magnetic mass mmag (inverse magnetic screening length) in 
the limit mmag ^ 7r by the following self-consistency condition: 

^' = ^^ 

mmag itself is given at high temperatures by 

"^mag = CN.g^T , (17) 

where c^^ is a perturbatively uncalculable number which so far can only be determined 
from lattice QCD calculations. 

The further evaluation of (15) is done differently in Refs. 35, 36. In both papers 
the approximation 7^ = 7^ = 7 is assumed, but in Ref. 35 Eq. (16) was approximated 
by the mmag independent expression 

7 = -a,iVer/2 Ina, , (18) 

while in Ref. 36 Eq. (16) was expanded in powers of 7/minag 

(19) 

where rj = corresponds to leading order and rj = iVc/(47rcjVe) to next to leading 
order. In the calculation c^^ = C3 = 0.405 was chosen^^. Further Altherr et al. 
used for Renr(g) = 3m^/2 the asymptotic limit for large q and go while Bilic et 
al. evaluated the dispersion relation (9) numerically in the high temperature (hard 
thermal loop) approximation. 

The results of both calculations are shown in Figure 3, with g = 2 and Nc = 3. 
Figure 3a shows a clear dominance of the thermal gluon decay in the mass region 
relevant for strangeness production [nis/T < 1). But the other two channels are prob- 
ably underestimated. They were calculated in the limit M ^ T from (3) and (4), 
respectively. An exact numerical treatment^^ of Eq. (3) and (4) shows a higher rate 
for these processes as seen in Figure 3b. The solid line corresponds to the sum of 
light quark-antiquark scattering and gluon fusion; the gluon decay is shown as dashed 
lines for various approximations for the damping rate. We see that the qq and gluon 
fusion processes still dominate the strangeness production rate. According to Ref. 36, 
the gluon fusion process is the more important one of the two channels, exactly as in 
the original tree-level calculation^^. Note the large influence of the so far unknown 
gluon damping rate on the gluon decay into ss pairs: a large damping rate, i.e. a 
large gluon width, enhances the strangeness production rate drastically. In view of 



In 



"^mae + 2m^ag7r, 



+ 1.1 



(16) 



7r = {1 + v) 



In 



+ 1.1 



mt 



mae / 
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Figure 3a: Strange quark production rates 
for different channels according to the calcu- 
lations of Ref. 35, from where this figure was 
taken. For details see text. 



Figure 3b: Strange quark production rates 
for thermal gluon decay with different damp- 
ing rates 7. The upper short dashed line 
corresponds to the damping rate of Ref. 35, 
Eq. (18), the next lower dashed line to 77 = 
and the lowest dashed one to 77 = Nd {'^t^c^^) 
in Eq. (19). The solid line is the sum of the 
gluon fusion and quark anti-quark collision 
rates. Figure from Ref. 36. 



this sensitivity the approximation 7^ = 7l should be reinvestigated and, if necessary, 
improved. 

The physically most relevant quantity is the total strangeness equilibration time 
Ts. The additional channel of thermal gluon decay can only shorten r^, but quantita- 
tively its influence is not really very drastic^^"^^. For temperatures around 200 MeV, 
Ts is of the order of > 10 fm/c; it strongly depends on T, the strange quark mass and 
the coupling constant a^. The controversy about the importance of the thermal gluon 
decay does not affect the conclusion that it won't be able to reduce the strangeness 
equilibration time to values well below 10 fm/c. Thus, even in a heavy-ion induced 
QGP, full strangeness equilibration remains questionable. Calculating the thermal 
gluon decay is a first attempt to include non-perturbative features in the strangeness 
production. At temperatures around 200 - 300 MeV such non-perturbative aspects 
play an important role, and therefore many more detailed investigations of this ques- 
tion (in particular the consistent implementation of the real and imaginary parts of 
the thermal propagators also in the diagrams la,b^^) are necessary to obtain a reliable 
picture. 
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2.2 Hadronic Production Mechanisms 

It is generally assumed that with the beam energies available at the AGS and SPS we 
may reach energy densities close to, but certainly not very far above the deconfinement 
phase transition. It is therefore a possibility that in the present heavy-ion experiments 
the collision zone spends all or a large fraction of its lifetime in the confined hadronic 
phase. With the critical temperature and energy density from lattice QCD now 
apparently having settled firmly42 Tc<150 MeV and ee<l GeV/fm^ we must, 
however, expect some deviations from the often used picture of a non-interacting 
mixture of hadronic resonances'*^. Being close to the point of deconfinement and 
chiral symmetry restoration, the masses and widths of the hadrons as well as their 
interaction cross sections may already be strongly modified by the medium compared 
to their vacuum values. 

Most authors working on this problem agree that the masses of most of the hadrons 
decrease with density and temperature'*'*"'*^. The Goldstone bosons of chiral sym- 
metry breaking (i.e. the pions and kaons) are an exception: their character protects 
their masses from large changes with temperature'*^''*^ and density^^. The reason for 
that is an approximate compensation of two medium effects. One can view these par- 
ticles as bound states of constituent quarks with a very large negative binding energy. 
This results in a very low total energy, i.e. mass of these particles. The constituent 
quark mass is generated by chiral symmetry breaking. By heating the system the 
constituent quark mass decreases, because chiral symmetry gets restored, but so does 
the absolute value of the binding energy. Altogether the Goldstone masses are nearly 
stable, and if anything their masses tend to slightly increase'*^ as the chiral phase 
transition is approached. 

Another controversial particle is the p meson whose mass some authors'*^'^"^ expect 
to drop in the medium while others^^ expect it to rise. Decreasing masses of heavier 
particles and resonances open the phase space for their production. This leads to the 
expectation that the mostly heavy strange particles are produced with larger rates in 
hot and dense hadronic matter^^ compared to calculations using free cross sections 
and mass thresholds. 

The effect of the medium on the hadronic production of strange particles has been 
studied mainly by Ko and his group; their work includes the production of kaons^'*, 
(f) mesons^^, and A's and A's^^. 

Kaon production 

Hadronic kaon production in heavy ion collisions proceeds mainly through the follow- 
ing three production channels [M: non-strange meson; B: non-strange baryon) 

(a) BB — > KYB 

(b) MB — > KY (20) 

(c) MM — > KK . 

The channel (a) is the dominant production mechanism in collisions of protons with 
hydrogen or light nuclear (e.g. Be) targets. There the if"'"/7r"'" -ratio at AGS energies 
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(10 - 14 A GeV incident momentum) is about^^ 0.04 . For heavier targets (e.g. Au) 
the ratio rises at the same incident momentum to-^^^ /ir~^ 0.1 , indicating that 
probably also channel (20b) begins to contribute. Channel (20c) will become impor- 
tant once a substantial pion density has emerged. This is supposed to be the case 
in heavy ion collisions, such as Si+Au, where indeed /ir~^ 0.2 at the AGS"^^^. 
Semi-hydrodynamical^ calculations^^ and also some Monte Carlo calculations^^, us- 
ing free strangeness production cross sections from experiment, fail to reproduce the 
if "'"/tt "'"-ratio of the AGS Si+Au experiment. The event generators for nuclear colli- 
sions either reproduce the pions correctly, underestimating somewhat the yield, 
as e.g ARC^°, or correctly reproduce the kaon yields at the expense of slightly over- 
predicting the pions, as RQMD^^. So far these calculations have ignored medium 
effects on the involved interactions; it is conceivable that somewhat larger elementary 
cross sections might give better agreement with the data. 

Ko et al. have parametrized the dependence of the hadronic masses on T and ps 
(baryon density) by^'* 

^* r /T'\2i"r 

(21) 

where the detailed behavior is model dependent. Tc is the deconfinement temperature 
(taken as Tc = 195 MeV in the studies described below), po is the nuclear ground state 
density [po = 0.17 fm~ ), and A is an empirical parameter, determined as A 0.2 - 
0.4 from K , p, and e~ scattering on nuclei^^. The power n depends on the model, in 
which the temperature dependence is studied, n = 1/6 is obtained from QCD sum 
rule calculations'*^, while n = 1/2 is favored by chiral models'*^. 

The density dependence of the kaon mass^^'^^ differs from the mass scaling (21). 
Within linear chiral perturbation theory the kaon mass behaves like^^ 

1/2 

(22) 



m* 


['-©1 


n 




(5)1 


m 











niK \ 

where pc = fK^n^K jTi^^ is the critical density for kaon condensation [fx = 93 MeV is 
the kaon decay constant, Ti^^ = 45 MeV is the kaon-nucleon sigma term). Eq. (22) 
reflects an exceptional behavior of the kaons: although the kaon belongs to the Gold- 
stone octet in a SU(3)/ representation for the flavor group, its mass decreases with 
density in contrast to the density scaling of Goldstone bosons discussed before. Ac- 
cording to the authors^^ the reason is the large explicit chiral symmetry breaking 
by the strange quark mass. This issue is, however, still under debate. Within the 
Nambu-Jona-Lasinio model the tendency for kaon condensation (22) is not supported: 
Lutz et al. get an increase in mass for the if"*", while the K~ stays nearly unchanged 
as a function of density^^. More recent calculations in the framework of chiral per- 
turbation theory lead to a rising and a decreasing K~ mass^'*; including higher 
orders of the approximation within a coupled channels approach reduces somewhat 
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Figure 4: Temperature and density de- 
pendence of {(Tv) for TTTT —f KK (dotted 
curve), 7rp KK (short-dashed curve), 
and pp —f KK (long-dashed curve). The 
solid curve is for MM —f KK, averag- 
ing over initial meson distributions. Note 
the different vertical scales in the upper 
and lower row. The figure was taken from 
Ref. 54. 



the decrease of the K~ mass , and additionally coupling in the A* (1405) resonance 
essentially destroys any indications for K condensation at high baryon densities^^. 

Thermally averaged cross sections (crv) of the strangeness producing reactions are 
needed in the calculation described below^'*. For a process 1 + 2 — )• 3...n they are 
defined by 

{o-i2^3...nVi2) = j d^ki j d^k2 fl{ki)f2{k2) ai2^3,„n{kl,k2)vi2 , (23) 

where fi denote the (normalized) thermal distribution functions for the incoming 
particles, and Vi2 is their relative velocity. The elementary medium cross section 
<^i2^3...n(^i5 ^2) is calculated using established form factors and coupling constants, 
but masses which satisfy the scaling relations (21) and (22). An example is shown 
in Figure 4 from Ref. 54, where the temperature and density dependence of mesonic 
production cross sections are plotted. One sees a clear rise of the cross sections with 
temperature and baryon density. 

The medium cross sections may now be used to calculate the final strange particle 
abundances in a heavy-ion collision. This requires a model for the time evolution of 
the nuclear collision. In Ref. 54 the hydrodynamical model of Biro et al.^^ was used. 
It employs an ideal hadron gas equation of state without phase transition. In addition 
to the hydrodynamical equations, which determine the overall time evolution of the 
initial fireball (carried mainly by pions and nucleons), rate equations for the evolution 

*The expression semi-hydrodynamical will be defined below. 
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Figure 5: The rapidity distributions of 
different particles at AGS. The figure is 
from Ref. 54 where the calculation (vari- 
ous curves) was made, while the data are 
from Ref. 128. 



of the strange particle densities were solved. The rate equation for a particle species 
s has the general form 



1 djVp,) 
V dt 



p1p]ps n Pi 
_ f 

^flnal-l 

^ pipjp° n pf ^ 



(24) 



The summation is meant to extend over all considered processes for producing (elim- 
inating) particles of kind s. The superscript indicates the equilibrium densities. 
Eq. (24) contains gain and loss terms and takes care of the expansion of the fire- 
ball volume V. Since hydrodynamical models for nuclear collisions usually assume 
complete chemical equilibrium (even for strange particles) throughout the dynamical 
evolution, we use the expression semi-hydrodynamical for these calculations which 
decouple some particle species from the bulk evolution and treat them separately. 

As initial conditions for Si+Au collisions at 14.6 A GeV a baryon density of 
Apo and an energy density of 1.4 GeV/fm^ were chosen, corresponding to an initial 
temperature of about 180 MeV. Figure 5 shows the rapidity distribution of pions, 
kaons and protons after freeze-out. Freeze-out is assumed to occur spontaneously at 
a common temperature of Tf = 125 MeV. 

The agreement with the data is reasonable; the kaon abundances are well repro- 
duced. The authors conclude (in contrast to the earlier work in Ref. 58) that no 
QGP phase transition is needed in order to understand the kaon yields at the AGS. 
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It should be stressed that the main mechanism for kaon production in the described 
analysis is the process pp — )• KK^^. The assumed reduction of the p-mass gives the 
necessary high initial p-density; however, as pointed out above, this behavior of the 
p must still be considered as controversial. 

(j) production 

The strangeness enhancement from p+A to A+A collisions is also clearly seen in the (j)- 
meson yield. The NA38 collaboration reported a rise of the ^/(/) + a;)-ratio from p+U 
to S+U by a factor of 3 in the dimyon spectrum (most central coUisions)^^. As origi- 
nally suggested by Shor^° it is possible to explain this enhancement as a consequence 
of QGP formation^^. However, hadronic rescattering mechanisms with different ab- 
sorption cross sections for the uj and (f) mesons by the nuclear medium^-^'^^ can 
also account for the observed systematics of NA38, although the required rather large 
fireball densities and lifetimes present some problems, on the conceptional level^^ as 
well as in practice, since the expected significant rise of the p dilepton yield with 
the fireball lifetime^^ is in contradiction with the observed independence of the pjuj 
dilepton peak with collision centrality^'*. 

There are also theoretical studies predicting a drastic mass reduction of the (j)- 
meson^^ in the medium. In combination with reduced kaon masses it is speculated 
that this could drastically enhance (j) production^^. The main hadronic production 
mechanisms are here 

(a) KK ct>N; 

(b) KK (f)p; 

these are the only channels taken into account in Ref. 55 in the rate equation (24) 
for the (j). The channel MM — )• KK violates the OZI (Okubo-Zweig-Izuka) rule and 
is therefore suppressed. 

The calculation of Ref. 55 takes similar initial conditions as in Sec. 2.2, but assumes 
a higher freeze-out temperature of Tf 140 MeV at CERN SPS energies. The result is 
a ratio of (f)/uj 20% which was claimed^^ to be in agreement with the measured one 
of NA38^^ in 200 A GeV S+U collisions. But a direct comparison with data obtained 
from the leptonic decay channel is problematic: First, the p and u> cannot be resolved 
experimentally; therefore the p should be included in the theoretical calculation to 
predict the (j)/{p + uj) ratio. Secondly, dilepton spectra are integrated over the whole 
history of the collision (which is important in particular for the dilepton yield under 
the p peak^^) while the calculation in Ref. 55 considers only ratios at freeze-out. Thus 
hadronic mechanisms for the observed (j) enhancement remain a subtle issue. 

It is expected that in the medium also the widths of the particles change. The 
calculation of Ref. 69 within the vector dominance model shows a large broadening of 
the ^- width with increasing density, and in addition (j) mesons and kaons have reduced 
effective masses. Contrary, the authors of Ref. 75 claim that within the Nambu-Jona- 
Lasinio model the width of the (j) decreases with increasing temperature, because 
the (j) mass goes down, while the kaon mass stays nearly unchanged, such that the 
decay phase space for the dominant channel (j) — )• KK closes. On the other hand, an 
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increase of the particle widths in the medium by collision broadening is a generic 
feature which should be taken into account in any realistic calculation. 

In the dilepton spectra no change in mass or width of the (j) has been seen so far 
within the (rather poor) experimental resolution. In the invariant mass spectra of 
K^K~ pairs from Si+Au collisions at the AGS an indication for a small mass shift 
of the (f) peak has recently been reported-^^^; the observed shift amounts to a mass 
decrease by 3-5 MeV for low-pr ^'^'s from very central nuclear collisions. A surprising 
feature of this reported mass shift is that in the selected kinematic region all ^'s 
appear to be shifted, i.e. the width of the shifted (j) peak is smaller than the mass 
shift and shows no traces of an unshifted component. Clearly this situation requires 
further clarification. 



A-A production 

In the context of strangeness enhancement as a QGP signal, the special role of 
strange antibaryons was pointed out very early-^-^. For producing antibaryons an 
energy threshold of at least 2 GeV must be overcome (taking free hadronic masses). 
The observed high A yield in 200 A GeV collisions is thus hard to explain by ordinary 
hadronic mechanisms. Hydrodynamical calculations^^ based on a hadron gas equa- 
tion of state with a freeze-out temperature*^ around 150 MeV fail to reproduce the 
A's. The same is true for event generators based on string models, unless string fusion 
(the formation of color ropes^^ or of quark-gluon clusters^^) is taken into account. 
Thus an explanation of the A yield on a microscopic level would give much insight to 
the nuclear reaction dynamics. 

Ko^^ et al. have investigated the A multiplicities from the NA35 experiment-^^^'-^^^. 
The analysis was done in the same manner as in Sections 2.2 and 2.2. An enhancement 
of A's is expected mainly due to a mass decrease with rising density and temperature. 

The mass scaling of the A (A) was deduced from the SU(3)/ Walecka model^° as 

= 1 2 Ps , (26) 

rriA mi. 



with the couplings (3/2) g^AA = QaNN = ^48 and a cr-meson mass of m^- = 550 MeV. 
Ps is the scalar nuclear density. 

The main contribution to the A production rate comes from meson collisions which 
are not suppressed by the OZI rule: 

(a) ^ AN; 

(b) Kp AN. ^ > 

As initial conditions'^ a temperature of To = 195 MeV and a density of twice the 
nuclear ground state density were chosen. The model results in a multiplicity of 1.5 A's 
per S+S collision (NA35, Ref. 157: 1.5 ±0.4 ; Ref. 165: 2.2 ±0.4), but underestimates 
the pion, kaon and A multiplicities by 15%, 30%, and 50% respectively. 



'Higher temperatures would solve the problem, but they are hard to justify. 
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We see that with certain initial condition and a mass reduction for the A (A) 
a reasonable value for the A multiplicity results. But the calculation does not give 
a consistent overall picture. Medium effects for the A's are certainly an important 
ingredient, but whether the mass reduction alone can account for the observed final 
yields remains to be proven. 

Further remarks on hadronic strangeness production 

The introduced model raises the question how the particles get on the mass-shell 
at freeze-out. At the chosen freeze-out temperature the particles still experience a 
considerable reduction of their free masses. Strictly speaking this is inconsistent with 
the naive freeze-out concept which assumes a more or less sudden transition from 
interacting to free-streaming particles. The refined picture behind these calculations 
assumes that after the "freeze-out" point all two-body collisions have ceased, but the 
particles still feel an attractive mean field which arises from virtual processes. It was 
argued^^'^-^ that on their way out the particles recover their vacuum masses by giving 
up kinetic energy as they move up the mean field potential well. In this process their 
momentum changes according to the equation 

m*^ + k*^ = ml + k\ (28) 

where the quantities with a star correspond to the in-medium values at freeze-out. 
This leads to an accumulation of particles with low asymptotic momenta. Shuryak 
has shown that such an assumption can explain the low-pj^ rise of pions^-^, and for 
the recently observed "cold" low momentum kaons at the AGS"^'*^ the same mech- 
anism was suggested^^. But there are also alternative explanations for the low-pj^ 
enhancement discussed in the literature^^, and in particular for the "cold" kaons the 
experimental situation is not yet completely settled-^^. Thus for the time being the 
direct experimental evidence for medium effects on the particles widths and masses 
remains shaky. 

The calculations of Ko et al. show that such effects are able to account for a large 
part of measured strange particle yields, without assuming QGP formation. In these 
calculations the chemical equilibration is driven by the strongly decreased masses, 
which was not taken into account in the earlier calculations of Ref. 33. However, 
up to now these investigations suffer from the fact that each calculation focuses on 
one particular aspect of the observations, without reaching a coherent and internally 
consistent description of all observables and their characteristic mutual relationships 
simultaneously. For example, an additional test for the model would be the multi- 
strange (anti-)baryons, measured by the NA36 and WA85 collaborations at CERN 
(for references see Table 9). In spite of the extremely high mass thresholds, their mul- 
tiplicities are high and present a particular challenge to any hadronic kinetic model. 

Other points need improvement, too: the study of Ko et al. suffers from a ther- 
modynamical inconsistency^'*. The reduction of the masses has only been taken into 
account in the particle distribution functions and cross sections for the kinetic rate 
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equations, but not in the equation of state for the hydrodynamic evolution. One ex- 
pects that the mass reduction should also lower the energy density and the pressure 
of the system. 

Last but not least the subject of medium effects on hadron widths and masses 
is altogether still very controversial. Different approaches lead to qualitatively dif- 
ferent and mutually contradictory results. A more fundamental and comprehensive 
approach to the problem of in-medium masses (based, for example, on QCD-inspired 
low-energy effective hadron Lagrangians) should eventually replace the rather phe- 
nomenological scaling law Eq. (21). Since these effective masses enter the chemical 
equilibration rates in a crucial way and the whole issue of strangeness equilibration as 
a signature for interesting new physics is mostly a question of time scales, this should 
be a high priority project. 

3 Thermal Models 

In this chapter we explain in more detail the physics of thermal models and their 
application to relativistic nuclear collisions, and discuss the conclusions that can be 
drawn by looking at strange particle production yields. 

3. 1 Thermal Parameters 

The thermal model starts with the assumption of a locally equilibrated, not neces- 
sarily stationary fireball. The fireball could consist of a QGP, a hadron gas, or an 
equilibrated mixture of both. The following set of thermal parameters can be used 
for all these phases. 

The fireball is described by a volume V , a temperature T and chemical potentials 
/i for the conserved quantum numbers. In general the intensive thermal parameters 
are local quantities, i.e. depend on space and time. This description corresponds 
to a grand canonical ansatz, i.e. the quantum number densities parametrized by the 
chemical potentials are not conserved exactly, but only in the average. The corrections 
compared to a canonical description with exact quantum number conservation are of 
order 1/ViV, with N being the number of particles with corresponding quantum 
numbers. Thus the corrections are expected to be negligible for heavy nuclei. 

On the time scale of a nuclear reaction the important forces are the strong and 
electro-magnetic interactions, both conserving the quark flavors. Flavor changing 
weak interactions can be completely neglected in nuclear collisions. (This is a first 
example of the entrance of time scales: in the cosmological context, where the dynam- 
ics happens on the time scale of microseconds, weak interactions can be considered as 
chemically equilibrated, and the corresponding parametrization in terms of chemical 
potentials and conserved quantum numbers is correspondingly different.) Therefore 
we introduce chemical potentials /iu, /id and /ig for the three flavors considered here. 
Note that this ansatz automatically takes care of baryon number and electric charge 
conservation. In view of the (approximate) isospin symmetry of the colliding nuclei, 
which carries over to the fireball, it is common to combine the two light flavors 

/^q = ^ (/^u + ^id) ■ (29) 
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The small breaking of the isospin symmetry can be parametrized by 

Sfi = fid - flu ■ (30) 

The nuclei considered here will be either exactly isospin symmetric (e.g. ^^S), or 
the breaking is small, leading to Sfi 0.^^ Therefore isospin effects can usually be 
neglected, and we use always the averaged light quark chemical potential (29). Flavor 
conservation is expressed on the level of hadrons as baryon number and strangeness (as 
well as isospin) conservation. It is also common to introduce the chemical potentials 
on this level, i.e. a baryon chemical potential fiB and a chemical potential for the 
quantum number strangeness fis- These two equivalent sets of chemical potentials 
are related by 

flq = /iB/3 , 

/is = /iB/3 - fis , 

/is = /iq - /is , (31) 

where the minus signs are due to the conventional assignment of strangeness -1 to 
the strange quark. 

If Sh are the baryon number and strangeness of hadron its chemical potential 
can thus be written either as 

fih = BhfiB + Sh fis (32) 

in terms of /tB and /ts, or as 

t^h = i^ht^q + yltis , (33) 

where z/^, z/^ count the number of light and strange valence quarks inside the hadron, 
respectively, with anti-quarks counted with a minus sign. 

The particle numbers are more directly given in terms of the fugacities, related to 
the chemical potentials by 

\i = e^'/^ . (34) 

The fugacity of each hadronic species is simply the product of the valence quark 
fugacities, viz. Ap = A^Aa, Ak+ = AuAg, etc. 

3.2 Relative and Absolute Chemical Equilibrium 

On the nuclear time scale, the state of absolute chemical equilibrium is defined as the 
state in which all inelastic processes which transform different hadron species into 
each other, subject only to the conservation laws for baryon number and strangeness, 
have come to equilibrium. On much larger time scales, this is actually only a state 
of relative chemical equilibrium, because weak flavor changing processes have not 
been equilibrated. Since the strong and weak interaction time scales are so widely 
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separated, the concept chemical equilibrium with respect to strong interaction pro- 
cesses only is very useful. In such a state the chemical potentials for particles and 
anti-particles are opposite to each other, implying for the respective fugacities 

k = K' ■ (35) 

However, the lifetime of a nuclear collision is actually so short that not all allowed 
chemical processes among the different species of hadrons may have time to fully 
equilibrate. The assumption of absolute chemical equilibrium is thus, even at freeze- 
out, the most questionable one in a thermal fireball approach. Thermal equilibration 
is much easier to justify since every type of collision (elastic and inelastic ones) in- 
volves exchange of momentum and thus leads to thermalization of the momentum 
spectrum. No thresholds or small cross sections inhibit thermal equilibration. Chem- 
ical processes, on the other hand, involve only certain inelastic channels with some- 
times small branching ratios (cross sections) and may be suppressed by large mass 
thresholds (in particular in the case of baryon-antibaryon or strange-antistrange pair 
production). 

Looking at the systematics of hadronic chemical reactions^^ it is, however, possible 
to subdivide them into two classes with generically well-separated time constants: The 
first class of processes involves the creation or destruction of only light quarks (or 
hadrons containing only light quarks); these processes have low mass thresholds and 
happen relatively fast (a few fm/c). Such processes include the exchange of existing 
strange quarks or antiquarks among different hadron species, e.g. in the process 
K N — > ttY, which is characterized by a strong absorption cross section of antikaons 
by nucleons. Much slower are processes involving the creation or destruction of one 
or more strange quark-antiquark pairs (several 10 fm/c). A quite analogous situation 
is found in a QGP environment: qq pairs are created from gluons with a much faster 
rate than ss pairs from gluons or light quark pairs. 

Given these different time scales it is thus sensible to assume that the state of 
absolute chemical equilibrium with respect to strong interactions is reached in stages: 
first a state of relative chemical equilibrium is reached in which gluons and light 
quarks and antiquarks are thermally and chemically equilibrated, but strange quarks 
and antiquarks are only thermally equilibrated, but haven't completely reached their 
equilibrium density. On the hadronic level this state corresponds to chemical equi- 
libration of all non-strange hadrons as well as a partial chemical equilibrium with 
respect to all strangeness exchanging processes, while at the same time the absolute 
amount of strange valence quarks and antiquarks remains below its saturation limit. 
In other words, all fast chemical processes are equilibrated, while the slow strange 
pair production processes are out of equilibrium. Only in a second, much later stage 
also these latter processes equilibrate and a state of full chemical equilibration occurs. 

The crucial question is how these various chemical time scales compare to the 
lifetime of a nuclear collision fireball. As far as we know from the discussions in 
Section 2, this hierarchy between light and strange quark producing processes is 
generic and true in both hadronic and QGP environments; the absolute time scales 
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for both types of processes are, however, much faster in the QGP than in the hadron 
gas, both because of the larger densities of scatterers and lower mass thresholds. This 
is the deeper reason why strangeness is such a useful observable: with an initial QGP 
phase, relative chemical equilibrium in the above sense should be a valid concept 
in any case, and the nuclear collision lifetimes may even be at the threshold for 
complete chemical equilibration. Without the QGP as a "catalyst" absolute chemical 
equilibrium seems out of reach in nuclear collisions, and even the concept of relative 
chemical equilibrium may not be very well realized. 

The idea of relative chemical equilibrium was first introduced heuristically by 
Rafelski^'* and later exploited in more detail by several groups (see Refs. 16, 25, 21, 
22, 24, and 17). Only recently it was put on a solid theoretical foundation and gener- 
alized to other situations (e.g. the possible breaking of chemical equilibrium between 
mesons and baryons after sudden hadronization of a QGP followed by rapid freeze- 
out^^) by C. Slotta^^. He used the principle of maximum entropy to derive the distri- 
bution function for relative chemical equilibrium. In the case of relative strangeness 
equilibration it turns out to be sufficient to introduce one common strangeness satu- 
ration factor 7s with < 7s < 1, where 7s = 1 means full chemical equilibrium resp. 
complete saturation of the strangeness phase space. The strange quark and antiquark 
densities are then regulated by the effective fugacities 

Af = 7sAs , 

Af = 7sA7i . (36) 

In the hadron language these enter the hadron fugacities according to Eqs. (33,34). 
The individual hadron fugacities Xh are thus given by a product of effective valence 
quark fugacities, 

A. = n Af , (37) 

j 

and the approach to equilibrium is controlled by the factor: 

7. = 7s"^''^^ , (38) 

where the power n[h, s) counts the total number of valence strange quarks plus anti- 
quarks in the hadron species h. From the treatment in Ref. 85 it follows that 7s must 
be treated as an additional fugacity, i.e. it can be written in the form 

7s = e^H/^ , (39) 

where /i|s| is an additional chemical potential for strangeness carrying quarks and 
equal for strange and anti-strange quarks. Its deviation from zero parametrizes the 
deviation of the thermodynamic state from absolute chemical equilibrium. The hidden 
temperature dependence on 7s in Eq. (39) is important for the calculation of the 
entropy of the system^^. 

Since 7s parametrizes the relative equilibrium between the two subgroups of non- 
strange and strange particles, it can be used as a quantitative measure for the 
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strangeness enhancement from p+p over p+A to A+A collisions . Due to its ther- 
mal definition, however, this requires an application of the thermal model even to 
low-multiplicity p+p collisions. The viability of such an approach has been discussed 
for more than 30 years (see Ref. 43 for a historical overview and Ref. 86 for a recent 
careful investigation of the pp- annihilation system at rest) and can up to now only 
be justified by its phenomenological success. 



3.3 The Partition Function 
Partition function for the QGP 

The calculation of all thermodynamical variables starts from the grand canonical 
partition function Z. It is related to the thermodynamical potential Q, through 

^ = -T\nZ. (40) 

This can then further be expressed by the pressure P = —Q,IV. For the QGP with 
two massless flavors (u and d) and a massive (mg) strange quark we have 

T 

37 1 

90 27r2^'i^<i 

TTlg 

Here B is the MIT-bag constant which is needed to simulate the background confine- 
ment pressure when we investigate the strangeness phase space diagram. For other 
thermal quantities it has no relevance. From the partition function one derives the 
net particle densities for light quarks pq, strange quarks ps, and the baryon density 
Pb as 

= V [-^) = V i^ATj = 1 ' ^''^ 

= I (P^ + P^) ■ (44) 



Partition function for the hadron gas 

The partition function for the hadron gas is given as a product of the one-particle 
partition functions of the different hadrons h: 

^«(r,y,/Xq,/X3,7s) = ]l^xp[ZH{T,V,tiH,^h)] . (45) 

h 
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The latter is given by 



Qh is the degeneracy factor, mh is the mass, is the chemical potential according to 
Eq. (32,33), and 7/1 is the strangeness suppression factor of hadron species h given by 
Eq. (38). The upper index "pt" stands for point-like, because so far we treated the 
hadrons as point-like particles. The real pressure involves an additional proper 
volume correction^^ and is given by 

= Trivia <*^' 

where e^' is the energy density calculated for point-like particles. 

In the analysis of measured strange particle yields the Boltzmann approximation is 
sufficient. Quantum statistical corrections are important only for the very light pions 
or for fermions at very large baryon densities (which will not occur in the actual 
applications discussed here). The emphasis of a chemical analysis of the data lies on 
the fugacities and the strangeness suppression factor 75. We spell these factors out 
explicitly and obtain in Boltzmann approximation 

In^" = ^ [Fm + (A^ + A-^)Fn,a + (AsA-^ + \:^\)-i,F^ (48) 
+(A3A^ + K'K'HF^ + (A^^A, + K'X-'KF^^ + (A^ + K'HFn] . 

This expression takes into account all non-strange mesons (M) (for mesons with hid- 
den strangeness see below), kaons (K), non-strange baryons, i.e. nucleons (N) and 
deltas (A), hyperons (Y), cascades (S) and the omegas (Q), together with their anti- 
particles and all their resonance excitations. The factors 

Ff = Y.9H{^\MmHlT) (49) 

h 

(where K2 is the modified Bessel function of second order) result from a summation 
over all hadrons (ground state and resonances) h in the hadron family / where / = 
M, N, A, K, Y, S, or Q. 

In order to describe consistently a strongly interacting hadron gas, the full (diver- 
gent) hadron spectrum should be taken into account, as described by Hagedorn'*^ . For 
a comparison with data we need, however, the final decays into stable (measurable) 
particles. These branching ratios are only known for the lowest states. Therefore one 
must cut the hadronic mass spectrum somewhere. The bias from this artificial cut 
will be investigated when we analyze the NA35 data from S+S collisions. 

The effect of the strangeness suppression factor 7s on mesons with hidden strange- 
ness via Eq. (38) requires a careful discussion. According to the SU(3)f quark model 
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classification^^ we take (/>, /o(980), /i(1510), /2(1525), (/)(1680) and (/)3(1850) as pure 
55-states, with a suppression factor 7^ . The physical rj- and 7/'-states are mixed states 
consisting of uu, dd and ss. We take a mixing angle between rj and rj' of —20°^^. 
This results in an 55-content of 32% in the rj and 68% in the rj'. It follows that the 
multiplicities of these two mesons incorporate the following strangeness saturation 
factors: 

7^ = 0.68 + 0.3273' , 7^, = 0.32 + 0.6873' • (50) 
Therefore the term Fu contains some 7s factors not explicitly written down in Eq. (48). 

4 The Phase Diagram of Strange Matter 

In this Section we describe the phase diagram for strange hadronic matter. For sys- 
tems with vanishing net strangeness it was first studied in Refs. 89, 90. Due to the 
conservation of strangeness by the strong interaction this is the most relevant case for 
heavy-ion collisions. Due to surface radiation of strange particles from the expanding 
collision fireball it is, however, possible that the net strangeness of the system fluc- 
tuates around zero or even acquires an appreciable non-zero value^^. Furthermore, 
hadronic systems with finite net strangeness are relevant for early cosmology and the 
structure of neutron stars (with or without a quark core) where weak interactions 
are equilibrated. Therefore the authors of Ref. 92 generalized the phase diagram to 
systems with non-zero net strangeness. 

Such systems can be characterized by their strangeness fraction 

_ ps{T, /Xq,/Xs) 

the net number of strange quarks (or minus the strangeness) per baryon. 

Strange hadronic matter, like matter with zero net strangeness, undergoes a phase 
transition to a quark gluon plasma if its density is sufficiently high. High particle 
density can be achieved either via quark pair creation by heating the system or via 
compression (high net baryon density). In the latter case it doesn't matter whether 
the baryon density is generated by strange or non-strange baryons, nor is the sign of 
the net baryon density important; only its modulus must be large. In the following 
subsection we discuss this phase transition. 

4.1 The Strange Matter Igloo 

Good model equations of state for hadronic matter near this phase transition are 
hard to obtain. The usual phenomenological procedure is to estimate the location 
of the phase transition by matching two different model equations of state, one for 
the hadronic phase which is valid at low temperatures and normal baryon densities 
where we have information from conventional nuclear physics, and one for the quark- 
gluon plasma above the phase transition which is derived from QCD at very high 
temperatures where it can be calculated perturbatively. Such a matching procedure 
produces almost by necessity a first order phase transition. 



23 



Two commonly used equations of state are in this context the ones defined in 
Eqs. (41) and (45)-(47), respectively. The phase transition line is then obtained from 
the Gibbs stability conditions for phase coexistence: 

Pb = Pq, Th = Tq , /iq,H = /iq,Q , and /is,H = /^s.q • (52) 

By setting the temperatures and chemical potentials in the two phases equal, one 
obtains the pressure balance relation 

PH{T,Hq,Hs) = Pq{T,Hq,Hs) . (53) 

Its solution defines a 2-dimensional phase coexistence surface in the 3-dimensional 
half-space defined by T (>0), fiq, and fis, the so-called "strange matter igloo"^°. 
Along the igloo surface hadronic and quark matter can coexist with variable volume 
fractions, but in general with different strangeness fractions and in the two 
subphases. If the hadronic subphase occupies the volume fraction a = Vh/(^ + Vq), 
the strangeness fraction of the total system depends on a and is given by 

fs{T,Hq,Hs) = af^{T,Hq,Hs) + (1 - a)f^{T,Hq,Hs) . (54) 

Systems with a fixed strangeness fraction fs will thus move along the igloo surface 
as they hadronize: as a increases from (QGP) to 1 (hadron gas), the condition of 
constant strangeness fraction (54) cuts out a strip from the igloo surface, as shown 
in Figure 6. If the system is further constrained to possess a certain fixed baryon 
number B and total entropy S and is allowed to hadronize slowly and adiabatically 
(i.e. at constant specific entropy S/B), it will follow a specific line, the so-called 
isentropic expansion trajectory, through that strip. 

4-2 Isentropic Expansion Trajectories 

A discussion of these trajectories is very interesting. Suppose that our system reaches 
a state of thermodynamic equilibrium somewhere in the quark phase (characterized 
by a point Pq in Figure 6 outside the igloo) and begins to expand hydrodynamically. 
In the absence of shock discontinuities the hydrodynamic equations conserve entropy. 
Since baryon number and strangeness are also conserved by the strong interactions, 
the system will expand at constant strangeness fraction fs and specific entropy S/ B. 
These two conditions together define a line through the three-dimensional phase dia- 
gram. This is the isentropic expansion trajectory just mentioned. It begins at Pq in 
the quark phase, reaches the igloo surface at a point Pi at the edge of the strip on 
that surface which corresponds to the chosen /s, and then begins to cross that strip. 
Along the trajectory both the temperature and baryon density are continuously de- 
creasing functions. If allowed to stay in local thermal equilibrium forever (in reality 
the system will decouple once the density becomes too low), the hydrodynamic ex- 
pansion will only stop at zero temperature and zero baryon density, after the system 
has expanded over an infinite volume. There is no other possibility to reach a state 
of zero temperature at finite total entropy. 

The interesting question is which point on the T = plane the trajectory will 
finally approach. There are two possibilities: 
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Figure 6: The phase diagram of strongly interacting matter in (T, /Xq,/is) space for B^l^ = 180 
MeV. The igloo-type surface describes the phase-coexistence region (mixed phase) between quark 
matter (outside) and hadronic matter (inside). Figures (a) through (f) show various sections 
through this surface corresponding to systems with fixed strangeness fraction /g. Also shown are 
the projections of the mixed phase region onto the three coordinate planes. From Ref. 92. 



1. The trajectory reaches the other side of the phase-coexistence strip before hit- 
ting the T = plane. In this case it will leave the igloo surface at a point P2 
on the other side of the strip and proceed into the interior of the igloo until it 
finally reaches a point P3 on the T = plane. 

2. The trajectory reaches the T = plane before the crossing is completed, at a 
point P3 somewhere on the line in which the strip and the zero temperature 
plane intersect. 
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These two cases correspond to two completely different physical situations. In 
order to discuss them we must understand the T = limit of the phase diagram in 
more detail. 

4.3 The T — )• Limit of the Phase Diagram 

Non-strange hadronic matter has at zero temperature a quite simple structure: mesons 
are absent (because they don't carry any conserved charges and thus their chemical 
potentials vanish), and only those baryons exist whose mass is smaller than the baryon 
chemical potential fiB (i.e. the Fermi energy). If fiB is smaller than the mass of the 
lightest baryon, the nucleon, then no particles exist at all. The hadronic state corre- 
sponding to T = and fiB < m^ is thus the vacuum; it has zero pressure. In other 
words: as T — )• 0, any state with finite net baryon number has to approach a point 
with fiB > mjv- 

If strangeness is added, things get more complicated. The T = phase diagram 
is now a 2-dimensional plane spanned by fi^ = /iB/3 and fis. Its structure (phase 
coexistence line, lines of constant strangeness fraction, etc.) is very rich and highly 
nontrivial. A detailed discussion is given in Ref. 92. Here we can only summarize the 
most important features. 

Generally speaking, mesons are still absent at T = 0. But there are exceptions: 
since kaons are bosons which carry a conserved charge (strangeness), they can form 
a Bose condensate if their chemical potential is equal to their rest mass: 

/^if = /iq - /^s = rriK or ^g- = jj,^ - jj,^ = mx ■ (55) 

These conditions define two straight lines in the /iq-/is plane, where the first corre- 
sponds to kaon, the second to antikaon condensation (Figure 7). These equations 
delineate the physical limits of the hadron phase at T = which can only exist 
between these two lines. 

Baryons exist at T = only if their chemical potentials exceed their rest masses. 
The baryon chemical potentials have the form fii = Sifis + (3 — Si)fiq where Si is 
the number of strange quarks in the baryon. For 5^=0, 1, 2, 3 (nucleons, hyperons, 
cascades and Q-baryons) the condition fii > mi defines four straight lines. Their 
envelope separates the region with zero baryon density (and zero pressure!) from the 
region with finite baryon density (see Figure 7). 

One can now solve the pressure balance equation (53) in this plane. For the three 
values 51/4 = 145, 195, 235 MeV for the QCD bag constant one obtains the three 
solid curves shown in Figure 7. One sees that for the two smaller bag constants 
part or all of the coexistence line lies inside the region of zero baryon density in the 
hadronic phase; these parts of the line thus correspond to a mixed phase, where the 
hadronic subphase is just the vacuum. Another way of describing such a state is 
that it consists of quark matter bubbles separated by vacuum! This state is stable, 
because it solves the pressure balance equation: the pressure of the hadronic phase 
(=vacuum) is zero, and the quark bubbles have also zero pressure, because the Pauli 
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Figure 7: Phase structure at T = 0. The 
three solid curves are the phase coexistence 
lines for B^/^ = 145, 195, 235 MeV. The 
thick solid and dashed lines separate the re- 
gions of zero baryon density (left) from re- 
gions of finite baryon density (right). The 
left region, which at T = is empty, is sub- 
divided into 6 smaller regions, indicating 
which particle species dominate at small, 
but non-zero temperature. The small num- 
bers along the lines refer to equation num- 
bers in Ref. 92 from where this figure was 
taken. 



pressure of the quarks is balanced by the bag pressure on the surface. If the quark 
matter in this state has fs ^ 0, these bubbles are called strangelets. 

One can show^^ that for B^^^ < 200 MeV some part of the phase coexistence line 
lies in the region of zero baryon density in the hadronic phase. Actually, practically 
all points corresponding to a non-zero strangeness fraction lie on this section of the 
phase coexistence line. For these values of the bag constant cold strange quark matter 
is thus mechanically stable and stops expanding; it cannot hadronize. 

4-4 Hadronization of Strange Quark Matter at T ^ 

Another way of plotting these results is shown in Figure 8. Here we show the phase 
diagram in the (T, pe) plane, eliminating fi^ and fis. The previous phase coexistence 
line is here resolved into a whole (shaded) domain, the mixed phase (M). Different 
points in this mixed phase at the same temperature T correspond to different volume 
fractions a of the hadron and quark subphases. 

One sees that, for the particular value of the bag constant chosen in this Figure, 
at low temperatures the mixed phase extends all the way to zero baryon density even 
for systems with very small amounts of net strangeness. This is just the situation 
discussed at the end of the previous subsection. If we start at finite temperature 
and let the system expand adiabatically, we see that some of the isentropic expan- 
sion trajectories never exit the mixed phase. This corresponds to the situation 2 
mentioned at the end of Section 4.2. Strangelets form. On the other hand, if the 
specific entropy S/ B was sufficiently large to begin with, this does not happen and 
the system hadronizes completely; this corresponds to situation 1 mentioned at the 
end of Section 4.2. 
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Figure 8: Isentropic expansion trajectories in the (T, pb) plane, for 5^/4 = 180 MeV and several 
/s values. For /s 0, at T = the mixed phase extends to zero baryon density. Systems with 
sufficiently large /s and/or small specific entropy SjB reach T = without leaving the mixed 
phase, forming strangelets. The full and open circles indicate pion freeze-out for spherical fireballs 
with 4 and 8 fm radius, respectively. For B^/'* > 200 MeV, the diagrams look like Figure (a) for 
fflZZ values of /s, and strangelet formation is impossible. From Ref. 92. 



It is, of course, interesting to ask which are the largest allowed values for SjB 
which still lead to strangelet formation in a system with given /g. We already 
know that we need bag constants below 5^/^ = 200 MeV {B < 200 MeV/fm^) 
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for strangelets to form at any value of SjB and f^. For 5^/^<145 MeV {B < 60 
MeV/fm^), on the other hand, even non-strange quark matter would be stable, which 
contradicts nuclear phenomenology. Within this narrow window of the bag constant, 
the values of fs and S/B play a crucial role^^. E.g. for B^^^ = 180 MeV, strange 
matter with /s=l (hyperon matter) completely hadronizes as soon as S/B > 6. Typ- 
ical values for S/B from heavy ion collisions at Brookhaven or CERN range from 13 
to 45 (see Section 7). 

We would like to close this section with one important remark: if (strange) 
quark matter undergoes adiabatic hadronization at fixed specific entropy S/B and 
strangeness fraction fs (which may be zero), the chemical potentials usually change 
during the phase transition. Thus a strangeness neutral QGP will start to hadronize 
with fis = 0, but at the end of the hadronization process fis will take on a differ- 
ent value which depends on T and fiB and corresponds to the equilibrium value in 
a strangeness neutral hadron gas. To draw conclusions from measured values of the 
chemical potentials in the final hadronic state about their values in the preceding 
QGP state is thus highly nontrivial. In particular, a vanishing strange quark chem- 
ical potential, fis = 0, in the final state can have survived from a primordial QGP 
phase only if hadronization did not happen adiabatically, but suddenly. 



5 Ingredients and Extensions of the Thermal Model 

5.1 The Strangeness Neutrality Condition 

An important constraint on the thermal model parameters arises from the fact that 
in strong interactions strange quarks are only created in pairs, and that therefore the 
nuclear collision fireball (up to small fluctuations due to possibly asymmetric surface 
radiation of strange and antistrange hadrons^ or due to "strangeness distillation" 
during hadronization of a QGP^^) always remains approximately strangeness neutral. 
The condition of vanishing total strangeness according to Eq. (43) takes the form 

= ^3 = (^)-(^) = A3^1n^. (56) 

For the hadron gas Eq. (56) is an implicit equation relating Ag to Aq in a way which 
strongly depends on the temperature^^'^^. 

In Figure 9 we show the relation between fi^ and fis for different temperatures 
which ensures vanishing net strangeness in a hadron gas. One recognizes a slight 
dependence on the number of hadronic resonances taken into account. The curves 
have a characteristic behavior and, for not to high temperatures, intersect the fis = 
line twice. 

At high temperature the heavy resonance states become important. Unfortu- 
nately, the high mass resonance states are poorly known, especially in the mesonic 
sector where the states above about 1.5 GeV become very broad and are difficult to 
identify experimentally. This triggers an artificial infiuence of the known resonance 
mass spectrum on the strangeness neutrality relation between mesons and baryons at 
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Figure 9: Strangeness neutrality curves for a hadron gas for different temperatures. The dashed 
line corresponds to T = 150 MeV, the solid line to T = 200 MeV, and the dotted line to T = 250 
MeV. In the left figure the hadron mass spectrum was cut at rrXcut = 1-6 GeV, in right figure at 
meut = 2.0 GeV. 



high temperature. We believe that the strangeness neutrality curves are unreliable 
above temperatures of about 250 MeV. 

Generally Eq. (56) relates the three parameters Ag, Aq and T (resp. /ig, /iq, and T). 
If two of them are given, the third can be determined from this equation which thus 
reduces the number of independent parameters. The strangeness neutrality condition 
depends also somewhat on 7s, but this influence is weak because it really enters only 
through the multistrange hadrons. They carry a factor 7" (n > 1) which survives in 
Eq. (56) after the leading linear 7s-dependence has been factored out. 

For a QGP strangeness neutrality is a matter of the strange quarks alone, in- 
dependent of the amount of light quarks, because the strange quarks are no longer 
bound together with other non-strange quarks into hadrons. Strangeness neutrality 
in the QGP is thus given by /ig = independent of fi^. This can be easily derived by 
inserting Eq. (41) into Eq. (56). 

As we will see later, the experiments at 200 A GeV suggest a value for fis very 
close to zero. While for a strangeness neutral hadron gas this solution is only possible 
for T < 200 MeV and a particular value of /iq, a strangeness neutral QGP leads to this 
solution naturally and independent of fi^ and T. This is the characteristic difference 
in the thermodynamics of a hadron gas and a QGP. 

A more detailed discussion of this issue will follow later. In order to be able to 
discuss violations of the strangeness neutrality condition we define 
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as a measure for the net strangeness. 



5.2 Transverse Flow and Transverse Momentum Spectra 

If the system is in local thermal equilibrium, then the observed momentum spectra 
should reflect the (average) temperature. However, since a thermalized source which 
is surrounded only by vacuum must necessarily begin to expand, the thermal motion 
is superimposed in the spectra by a dynamical component arising from the collec- 
tive expansion. This is most clearly seen in the measured longitudinal momentum 
or rapidity spectra of hadrons from heavy ion collisions which for both Brookhaven 
and CERN energies are much broader than a thermal distribution from a stationary 
fireball^^'^'*'^^; the strong collective flow component here is partially due to incomplete 
stopping of the colliding nuclei and partially due to additional hydrodynamical expan- 
sion generated in the later stages of the collision. In the transverse momentum spectra, 
which look approximately exponential like a thermal distribution, the identification of 
a collective flow component is more involved; the observed convex shape of the spectra 
at low momenta, which does not agree with the concave behavior of a thermal distri- 
bution near pj_ = 0, can be associated both with transverse collective flow^^ and with 
non-thermal contributions from resonance decays after freeze-out^^. This ambiguity 
was systematically studied and resolved in Ref. 96; theoretical arguments based on 
the consistency of the observations with the freeze-out criterium^^'^^ and direct ex- 
perimental reconstruction of the contribution of A decays to the pion spectrum^^'^'*^ 
confirm the existence of a strong collective flow component ((/3f) — 0.3 — 0.4) also in 
the transverse momentum spectra. 

For a purely thermal distribution the m± spectra and the double differential dis- 
tribution dN/[dydm\) are given in the Boltzmann approximation by 

- ^^TO , cosh 7/p-'^(™-L<^°='^J'-^'') • (<^R) 
dml - 47r2' m^Ki(^^J ^ ^^^e ^m^l e , (59) 

where Ki(a;) is a modified Bessel function. If one fits the measured m_\_ spectra 
with the expression (59), an "apparent temperature" can directly be read off from 
the exponential fall off. This is a very common procedure which is easily applied 
by both theorists and experimentalists. But due to the problem of resonance decay 
and flow contamination just mentioned, care should be taken not to denote these 
slope parameters indiscriminately by "temperatures". At AGS and SPS energies 
the resulting apparent temperatures are always of the order of 200 MeV or higher 
which is well above Hagedorn's limiting temperature'*^ for a hadron gas and above 
the deconfinement temperature from lattice QCD'*^. This is true even though the 
contributions from resonance decays tend to steepen the spectra and thus reduce 
the apparent temperatures, as will be discussed in Section 5.3. It appears to be 
theoretically inconsistent to associate such a large temperature with the observed 
final hadronic state from heavy ion collisions. 
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Transverse collective flow provides a simple resolution for this dilemma. It is 
a necessary dynamical consequence of early thermalization of the collision zone^^. 
The amount of hydrodynamically generated transverse flow depends, among other 
things, on the initial energy (entropy) density and the equation of state. These are 
poorly known so far, and the resulting theoretical freedom must be removed by the 
experiment. 

Transverse momentum spectra from a collectively expanding thermalized source 
can be parametrized by^^ 



±,.,.A,/r*m,K, I. ^ii^^ , (60) 



where Io(a;) is also a modified Bessel function, i2f the freeze-out radius and L is some 
normalization length. p{r) stands for the transverse flow rapidity 

p(r) = tanh-^/3^(r) , (61) 

and the transverse flow velocity profile /3_l(t") is usually assumed to have the form 

/340 = a(;^)". (62) 

If the system expands longitudinally in a boost-invariant way^^, the double differential 
spectrum is simply given by 

dydm\ dm\ 

Eq. (60) represents an approximation which assumes a particular shape of the freeze- 
out hypersurface^^'^^. A detailed study within a global hydrodynamical evolution 
model^^'^^ shows that (60) represents well the transverse momentum spectra even if 
they are calculated with a more sophisticated kinematical freeze-out criterium. For 
not too heavy particles and not too low values of p_\_ the fit is also rather insensitive to 
the power a in Eq. (62)^^'^^'^^. Thus we take here for simplicity a constant transverse 
flow velocity, a = 0. 

For large m_\_ the spectrum (60) falls off exponentially, but with a slope which does 
not reflect the true temperature T of the locally equilibrated system at freeze-out, but 
rather an apparent "blue shifted" temperature Tapp which in the above approximation 
is simply given by^^'^^ 

In the case p — )• (/3f — > 0) Eq. (60) reduces the purely thermal spectrum Eq. (59). 
But since the experimental spectra fix only the asymptotic inverse slope (64) accu- 
rately, while at lower values of pj^ the details of the spectra (in particular the observed 
strongly concave shape of the pion spectrum) are governed by an intricate interplay 
of transverse flow and resonance decays, the true temperature and transverse flow 
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Figure 10: All fit pairs (T, /3f) compati- 
ble with the measured mx-spectra. Every 
point on these curves results in a good fit of 
the computed m± spectrum (including res- 
onance decays and transverse flow) to the 
respective measured^*'^^^ particle spectra. 
From Ref. 96. 



velocity cannot be easily separated by the fit procedure and must be determined by 
additional theoretical input. Figure 10 shows this ambiguity of the fitting procedure 
for the transverse momentum spectra of five different hadron species measured in 200 
A GeV S+S collisions by the NA35 coUaboration^^'-^^^. The experimental m± spectra 
can be reproduced within a large range of temperatures (100 - 200 MeV) by adjusting 
/3f accordingly. Hydrodynamical studies^^ show that dynamical consistency between 
the longitudinal and transverse momentum spectra and with the freeze-out kinetics 
restrict the allowed freeze-out parameters to a rather narrow domain around T = 150 
MeV and (/3f) = 0.3c. In detail there remains, however, certainly some model de- 
pendence to these statements, mostly related to the equation of state used for the 
hydrodynamical expansion. In spite of our theoretical prejudice towards freeze-out 
temperatures in the 150 MeV region^^''*^ we will therefore treat for the chemical 
analysis T as a free parameter. 



5.3 Resonance Decays 

At temperatures above 150 MeV hadronic resonance production becomes very impor- 
tant. After freeze-out the resonances decay into stable particles which are indistin- 
guishable from their directly emitted brothers and sisters. The measured distribution 
of a particle species i is thus given by 



dy dm\ dy dm\ ^ dy dm\ 
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with the direct thermal contribution given by Eq. (58) or (60). The contribution from 
decays of resonances i2— )-z + 2 + ... + n (with branching ratio bR^i into the observed 
channel i) is calculated according to^^ 

dN- r ft 
^ l,s,M jiYjml (66) 



M ( dN, 



R 



PIpI - [ME* - M^m^ cosh(y - y)^ \dY dMlJ 

Capital letters indicate variables associated with the resonance i2, ^/s is the invariant 
mass of the unobserved decay products 2, . . . ,n, and the kinematic limits are given 

by 

s_ = [Y^mk] , s+ = {M-miY, (67) 

\fe=2 / 

y± = yisinh"' — , (68) 



_ E*m^ cosh(y -y)± p^^p*' - mj smh'jY - y) 
^ " misinh^(y-y) + m2 ' ^^^^ 

E* = ^(M' + m^-s), p* = ^E*' - ml . (70) 

In Eq. (66) gn{s) is the decay phase space for the n-body decay i2— )-z + 2 + ... + n; for 
the dominant 2-body decay (assuming isotropic decay in the resonance rest frame) it 
is given by 

^^^'^ = d^^^'-"'^^- ^^^^ 

Forn > 2 see Ref. 83. 

The resonance spectrum dNR/dy dM\ entering the r. h. s. of Eq. (66) is in general 
itself a sum of a thermal or flow spectrum like Eq. (58) or Eq. (60) and of decay 
spectra from still higher lying resonances; for example, most decays of high-lying 
nucleon resonances into protons or anti-protons proceed through the A(1232) or its 
anti-particle {e.g. iV(1440) — > A(1232) + tt — > p + 27r). This has been taken into 
account in the calculations. Extracting the degeneracy factors and fugacities of the 
decaying resonances, we write shortly 

- t P'^ dN^ 

iVf = -iR\RNP = jrXr / dy / dml E 9r bR^i-^—^ , (72) 

^y-window ^m^"' ^ dy dm"^ 

with the sum now including only resonances with identical quantum numbers; if these 
quantum numbers agree with those of species z, the sum is meant to also include the 
thermal contribution. We will use this notation below. 

Between particles and anti-particles we have the relation 

Nf = ^r\-^' = . (73) 
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5-4 Particle Ratios 

A quantitative comparison of the experimental (strange) particle yields or ratios with 
the thermal model prediction must take into account the kinematic acceptance win- 
dows of the experiments. The measured ratio of two different particles a and b is in 
general given by 



Na -R -^dydm]_ 

Nb JJ dydm_i 



dN, 



(74) 



„ ■ dydrn^, 



Here Tli denotes the experimental kinematic acceptance window for particle species i. 
The advantage of considering particle ratios is the cancellation of the (poorly known) 
fireball volume Vat freeze-out (assuming that both particle species freeze out simul- 
taneously). Of course, in general one must consider the (different) resonance decay 
contributions to the numerator and denominator in (74) according to Eq. (65). A 
first, naive estimate can, however, be obtained by neglecting the resonance contribu- 
tions and further taking for dNi/[dydm\) the purely thermal distribution (58). Then 
the integration kernels of (74) are independent of the particle rest masses niafi- If we 
then choose 71a = Tib, i.e. cut the data to the same rapidity and m± window, the 
integrals over the Boltzmann factors cancel in the ratio^'*, which then simply reduces 
to a quotient of fugacities and 7s factors: 



Na 



est. 



^ . (75) 



With proper flow distributions (60) the integral kernels in (74) are no longer inde- 
pendent of the rest masses of the particles, and the simple estimate (75) must be 
corrected by a proper numerical calculation. Similar comments generally apply to 
the contributions from resonance decays. 



5.5 The Specific Entropy 

A useful quantity for studying the evolution of a nuclear reaction is the specific entropy 
S/ as already discussed in Section 4. If we knew S/ B and the temperature, we 
would be able to decide whether we are in the QGP phase or in the hadron phase, 
because of the characteristic high specific entropy of the plasma compared to the 
hadron gas at the same temperature. Thus it would be good to have a measure for 
S/ B. A rough estimate is provided by the following considerations. 

In the ultrarelativistic limit (T ^ m) the entropy per particle is independent of 
temperature and approximately S/N 4 (a little less for bosons, a little more for 
fermions). If the three charge states (+, 0, -), in which the stable final state hadrons 
arise, are equally distributed, we have iV'°' = 3/2(iV"'" + N~). The net charged 
multiplicity N'^ — N~ is given by the number of incoming protons. Therefore the 
total baryon number is approximately B 2(iV"'" — N~). Putting all this together 
one can write 

S S iV'°' l{N+ + N-) N+ + N- 

— ~ A . 21 1 — 3 ! f76) 

B B ~ |(iV+-iV-) N+-N-' ^ ' 
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Figure 11: The product Dq ■ S/B as a 
function of Aq for different values of Ag. 
Solid line: As = 1.00; dashed line: As = 
0.95; dash-dotted line: As = 1.05. The ad- 
ditional parameters are = 1, and T is 
fixed by the strangeness neutrality relation 
(56). The upper curves are for Dq (with- 
out resonance decays), the lower curves in- 
clude the decay of resonances after freeze- 
out. Figure from Ref. 25. 



Defining the specific net charge Dq as 



Dc 



N+ -N- 
N+ + N- 



(77) 



we can rewrite this as 

(S/B) -Dq^K ^Z. (78) 

Dq is an easily measurable quantity and most precisely accessible in emulsion exper- 
iments. 

Of course, in a realistic hadron resonance gas at T ~ 150 — 200 MeV the ultra- 
relativistic limit S/N ?a 4 is not realized. We can test the accuracy of the estimate 
Eq. (78) numerically. We introduce the notation Dq for the specific net charge in the 
hadron gas, before post-freeze-out resonance decays. In the upper curves of Figure 
11 we plot S/B • Dq for a hadron resonance gas at different values of Ag near Ag = 1, 
as a function of Aq. The temperature varies with Aq according to the strangeness 
neutrality relation (56) and actually covers a wide range over the Aq interval shown. 
Obviously the product S/B • Dq is indeed roughly constant as suggested by the es- 
timate (78), but the constant is approximately 5 rather than 3. The reason is the 
larger specific entropy of the heavier particles compared to the assumed value of 4. 
However, after adding the resonance decay contributions to the total charged par- 
ticle multiplicity in the denominator of (77), the lower set of curves in Figure 11 is 
obtained. It is surprising to see how the decays bring the approximate constant K 
again down to a value around 3 as in the simple estimate Eq. (78)! We thus have 
with Dq^ a reasonably accurate measure for S/B from a thermalized and chemically 
equilibrated hadron gas. For hadronic systems far away from chemical equilibrium, 
however, the estimate (78) should not be applied without further scrutiny. 
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6 Chemical Parameters from Strange Particle Ratios 

In this part we want to confront the thermal description with experimental data. 
This will be done in several steps: first, the input parameters for the thermal model, 
T, /is, /iq, 7s, (V) will be determined from a sufficiently large subset of the measured 
quantities. Due to the non-linearity of the equations to be solved (see below), it is 
not a priori obvious how large such a minimal set of observables must be. Second, the 
validity of the thermal ansatz will be subjected to a test by making predictions for 
additional observables for which data exist or can be measured in future experiments. 
Third, systematic deviations of the model predictions with the experimental data will 
be analyzed and interpreted. 

Before proceeding further we should comment on the determination of the freeze- 
out temperature Tf . It can be extracted either from the hadronic mj^-spectra^^'^^ or 
from the particle ratios-^^'^^'-^^. However, it should be stressed that these two proce- 
dures lead conceptually to different freeze-out temperatures: The momentum spectra 
determine (modulo possible ambiguities from collective transverse flow) the thermal 
freeze-out temperature, at which all collisions stop and the momentum distributions 
are frozen in. The temperature extraction from particle ratios relies on the assumption 
of chemical equilibrium and thus determines the point of chemical freeze-out where 
the number-changing inelastic collisions cease and the particle abundances are frozen 
in. Generically chemical freeze-out occurs prior to thermal freeze-out, but in excep- 
tional situations like the sudden hadronization of a QGP into a dilute, free-streaming 
hadronic system both freeze-out processes can occur simultaneously. 

We would like to stress that in the context of the following analysis temperature 
values far above Tc ~ 150 MeV generally pose a threat to its theoretical consistency: 
taking such temperatures seriously one is necessarily led to doubt the lattice QCD 
evidence'*^ for a deconfining phase transition around rc=150 MeV, and one must face 
serious questions about the applicability of the non-interacting hadron resonance gas 
picture underlying the analysis (at r=200 MeV a hadron resonance gas is already 
beyond its dense packing limit). 

6.1 Thermal and Chemical Parameters from the WA85 Experiment 

The WA85 experiment provides data on multistrange baryon and antibaryon produc- 
tion on sulphur-tungsten collisions at 200 A GeV (see Table 9). The WA85 collabora- 
tion has measured A, A, S~, S"'", and Q"'"; very recently (see WA85 contributions 
to Refs. 12, 13 and Ref. 188) they also extracted the spectra and yields of neutral and 
charged kaons. The data on strange baryons have been analyzed independently by 
two different groups, Refs. 25, 100 and Refs. 18, 19, with mutually consistent results. 
We will follow here the major steps of the work of Ref. 25 and comment on the other 
work later. 

We concentrate here on the dataset published in Refs. 185, 186. There we have 
the four ratios 

Re = ^ = 0.41 ± 0.05 for y G (2.3, 2.8) and > 1.9 GeV; (79) 
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Ra = J = 0.20 ± 0.01 for y G (2.3, 2.8) and > 1.9 GeV; (80) 

Re,a = ^ = 0.19 ± 0.01 for y G (2.3, 2.8) and > 1.9 GeV; (81) 

R2,A = ^ = 0-40 ± 0.04 for y G (2.3,2.8) and > 1.9 GeV. (82) 

Note that these ratios are all measured in the same kinematic region, and that only 
three of the four ratios are independent because of the relation 

-Ra Re a ~ -^H Re, A • (^3) 

First estimates for WA85 

The common kinematical region of the ratios in Eqs. (79)-(82) allows a first estimate^'* 
via Eq. (75). Specifically we obtain 

j^x^xi i!K'K^ . X 

This system of equations is easily solved for Aq, Ag and 

a: = ReR-a'-. (85) 

A^ = RaRe^] (86) 

73' = Re,aRe:a. (87) 

from which one obtains 

Aq = 1.47 ±0.05 ; A3 = 1.03 ± 0.05 ; 73 = 0.28 ± 0.03 . (88) 

As we will see in the following subsection, the first two values are roughly correct, 
but the estimate for 73 is wrong by about a factor 2 due to uncorrected resonance 
decays. 

Resonance gas analysis for WA85 

Incorporating resonance decays and allowing for collective flow needs a numerical 
treatment. Because for the particle ratios we need the chemical freeze- out temperature, 
which is not necessarily the same as the thermal freeze-out temperature reflected in 
the momentum spectra, we take Tf as a free parameter and investigate three cases: 

A: a thermal model without flow, /3f = 0, where the freeze-out temperature Tf is 
assumed to correspond directly to the value Tapp = 232 MeV deduced from the 
slope of the transverse mass spectra of high-mj^ strange (anti-)baryons"^^^'"^^^; 
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B: a model with a freeze-out temperature of Tf = 150 MeV, i.e. a value consistent 
with the kinetic freeze-out criterium developed in Ref. 97. This temperature 
also agrees with the QCD phase transition temperature from lattice Monte 
Carlo simulations. It entails a flow velocity at freeze-out of f3{ = 0.41; 

C: in order to maintain zero net strangeness in a hadron gas fireball at relative 
chemical equilibrium between strange meson and baryon abundances, we con- 
sider also freeze-out at Tf = 190 MeV with /3f = 0.20. 



The system of equations to be solved may be written as 



Rs — — ^ 
A 



m_|_ ^m^. 



m_|_ ^m*^. 



m_|_ ^m^. 



(89) 
(90) 
(91) 



where the iVf are only functions of the rest masses and the given temperature T, as 
can be seen from (72). The result of the numerical solution of (89)-(91) is given in 
Table 1. 

Since the probability to excite the various resonances contributing to the observed 
final particle abundances depends on Tf, it is not surprising that slightly different 
values of the chemical parameters result from the same data in the three different 
scenarios. Comparing the numerical results of Table 1 with the simple analytical 
approximation of Eq. (88), which contained no contributions from higher resonances, 
we can draw the following conclusions: 

1. The extraction of Aq and Ag according to Eq. (84) is only very weakly affected by 
resonance decays and by the origin of the slope of the mj^-spectrum (thermal or 
flow). (The absolute value of the associated chemical potentials does, of course, 
depend on the choice of the freeze-out temperature.) The reason is that in 
Eq. (89) (and respectively (90)) the first (respectively second term) in the sums 
occurring in the numerator and denominator completely dominates. Neglecting 
then the other terms leads to the same equation as in (84). Therefore the 
determination of Aq and Ag via Eqs. (85) and (86) leads already to a precise 
result. 

2. Since the S~/A ratio is the crucial ingredient for the determination of the 
strangeness saturation factor 7s, the effects from resonance decays and flow can 
be clearly seen in this number. From Eq. (88) we had extracted 7s = 0.28 it 0.03. 
This could be drastically improved by the simple observation that the S°, which 
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Table 1: Thermal fireball parameters extracted from the WA85 data ' on strange baryon and 
antibaryon production, for three different interpretations of the measured mx-slope. Resonance 
decays were included. For details see text. 





A 

"thermal" 
T = 232 MeV 
A = 


B 

"thermal & flow" 
T = 150 MeV 
/3f = 0.41 


c 

strangeness balance 
T = 190 MeV 
/3f = 0.20 


As 
fis/T 
fis (MeV) 


1.03 ± 0.05 
0.03 ± 0.05 

7 ± 11 


1.03 ± 0.05 
0.03 ± 0.05 

4 ± 7 


1.03 ± 0.05 
0.03 ± 0.05 

6 ± 9 


K 

fiB (MeV) 


1.49 ± 0.05 
0.40 ± 0.04 

278 ± 23 


1.48 ± 0.05 
0.39 ± 0.04 

176 ± 15 


1.48 ± 0.05 
0.39 ± 0.04 

223 ± 19 


7s 


0.69 ± 0.06 


0.79 ± 0.06 


0.68 ± 0.06 


e 


-0.22 


0.37 





S/B 


18.5 ± 1.5 
0.135 ± 0.01 


48 ± 5 
0.08 ± 0.01 


26 ± 2.5 
0.12 ± 0.01 



is approximately as abundant as the A, decays into A + 7, and this fast electro- 
magnetic decay cannot be resolved experimentally. Therefore about 50% of the 
observed A's result from S° decays which improves the naive estimate (88) by a 
factor of 2 to a value of 0.56. Taking into account the full resonance spectrum 
raises the value of 7s even further to values around 0.7 - 0.8, surprisingly close 
to full equilibrium, i.e. 7s = 1. 

3. We see from Eq. (88) and Table 1 that fis is small and compatible with zero. 
This confirms the results of Ref. 84 based on the simple estimates (85), (86) 
and implies that the WA85 data on strange baryon and antibaryon production 
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from 200 GeVA S+W collisions establish a vanishing strange quark chemical 
potential. This result is stable against large variations in the temperature and 
transverse flow velocities and thus insensitive to the production and decay of 
high mass resonances. The reason is that these parameters are extracted from 
particle/antiparticle ratios in which the resonance decays largely cancel. A 
vanishing strange quark chemical potential is natural in a strangeness neutral 
QGP. Note also that also the light quark fugacity is practically unaffected by 
the complex pattern of resonance formation and decay as well as by the amount 
of transverse flow in the spectra. 

4. Assuming that at freeze-out the hadronic system is in a state of thermal and 
relative chemical equilibrium and thus fully characterized by the parameters 
given in Table 1, all other hadron abundances can now be predicted and the 
strangeness neutrality can be checked. The second to last row in Table 1 shows 
that under these assumptions strangeness neutrality is only established in sce- 
nario C. Although it is not necessary that strangeness neutrality is satisfied 
locally - the conservation laws ensure strangeness neutrality only for the total 
fireball - it is hard to see how local variations from strangeness neutrality above 
|e| < 0.1 could be generated dynamically or by fluctuations. Thus we get into 
serious problems with scenarios A and B. 

One should keep in mind that the strange mesons play a crucial role in the 
strangeness neutrality relation (56); in order to resurrect, for example, scenario 
B (because it features a reasonably low freeze-out temperature), it would be nec- 
essary to break chemical equilibrium at least between the mesons and baryons. 
This would require the introduction of at least one further non-equilibrium pa- 
rameter (in addition to the global strangeness undersaturation factor 7s). As 
we will see below and in Section 7, entropy considerations point in a similar 
direction. Under these conditions the applicability of the concept of relative 
chemical equilibrium in a thermalized hadron gas state to the observed final 
state in heavy ion collisions at CERN energies appears highly questionable. We 
will, however, continue a little further down this line by asking whether the 
concept works at least within the subspace of baryons and antibaryons. The 
following subsection will show that, strikingly, the predictions arising from this 
picture are in qualitative agreement with experiment. 

5. A similar extraction of thermal and chemical parameters from the WA85 exper- 
iment was done in Refs. 17, 18, 19. Because Refs. 17, 18 deal with older data of 
the WA85 experiment we only want to comment on Ref. 19 which uses the same 
input as the calculation described above. The strategy in Ref. 19 is somewhat 
different. They extrapolate the ratios of Ref. 186 to pj^ = and compare these 
values to fully integrated thermal particle ratios. They also force the system to 
be strangeness neutral, using the strangeness balance equation (56) to uniquely 
determine the chemical freeze-out temperature. Their results are T = 190 it 10 
MeV, fiB = 240 lb 40 MeV and 7s = 0.7, in agreement with the scenario C 
presented in Table 1 above. 
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Table 2: Input (*) and predicted high-mx particle ratios near central rapidity (2.5 < y < 3.0) 
for 200 A GeV S+W collisions. The parameters for the three "thermal", "thermal & flow" and 
"balanced strangeness" scenarios are as specified in Table 1. 



NilNj 


A 


B 


C 






m]^' = 1.9 GeV 


"thermal" 


"thermal k flow" 


"strangeness 








balance" 




0.41 ± 0.05* 


0.41 ± 0.05* 


0.41 ± 0.05* 


A/A 


0.20 ± o.or 


0.20 ± 0.01* 


0.20 ± 0.01* 


E-/A 


0.19 ± O.or 


0.19 ± 0.01* 


0.19 ± 0.01* 


S+/A 


0.40 ± 0.04* 


0.40 ± 0.04* 


0.40 ± 0.04* 


A/p 


0.60 ± 0.06 


0.62 ± 0.06 


0.59 ± 0.06 


A/p 


1.2 ± 0.1 


1.2 ± 0.1 


1.2 ± 0.1 




0.53 ± 0.05 


0.29 ± 0.03 


0.45 ± 0.04 




1.1 ± 0.1 


0.60 ± 0.06 


0.94 ± 0.09 




0.10 ± 0.02 


0.05 ± 0.01 


0.08 ± 0.02 


n+/A 


0.42 ± 0.08 


0.23 ± 0.05 


0.36 ± 0.07 




0.85 ± 0.25 


0.85 ± 0.25 


0.85 ± 0.25 


p/p 


0.10 ± 0.02 


0.10 ± 0.02 


0.10 ± 0.02 


K-/K+ 


0.53 ± 0.08 


0.51 ± 0.08 


0.53 ± 0.08 




0.7 ± 0.2 


0.4 ± 0.1 


0.6 ± 0.15 


1.0 <p±< 2.5 GeV 








K°/A 


0.41 ± 0.08 


1.0 ± 0.2 


0.56 ± 0.11 


K°/A 


2.0 ± 0.4 


5.1 ± 1.0 


2.8 ± 0.6 



Predictions resulting from the WA85 analysis 

It is evident that the extraction of three thermodynamic parameters (Aq, Ag, and 7s) 
from three independent particle ratios cannot be considered a test of the generalized 
thermal model considered here, in particular since the extracted values depend on 
how the temperature is fixed in the analysis. However, with the thermal parameters 
now fixed according to scenarios A, B, or C, further particle ratios can be predicted 
and are given in Table 2. 

Some of these ratios can be compared with new data from 200 A GeV S+W col- 
lisions. WA85 has published the following (partially preliminary) values^^^'-^^^'-^^^'-^^°: 

Q+ 

Rq = -pr: = 0.57 ± 0.41 for y G (2.5,3.0) and > 2.3 GeV 
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Req = = 1-7 ± 0.9 for y G (2.5,3.0) and > 2.3 GeV 

Rka = ^ = 1.4 ± 0.1 for y G (2.5,3.0) and 1.0 < p± < 2.5 GeV 

Rka = = 6.4 ± 0.4 for y G (2.5,3.0) and 1.0 < p± < 2.5 GeV 

Rk = ^ = 0.60 ± 0.05 for y G (2.3, 3.0) and > 0.9 GeV 

The Q/Q-ratio agrees within error bars with the prediction, but the large experimental 
uncertainty limits its usefulness. In the model Q/Q is given by and thus should 
be 1 if the strange quark chemical potential vanished exactly. - i2so disagrees with 
the prediction: only for scenario A with the highest freeze-out temperature the two 
values overlap within their error bars. For larger 7s the situation would improve. 
However, again the experimental error bar is very large, and a definite conclusion 
requires better data statistics. - The K~ /K'^ ratio is much more accurately known; 
its central value is slightly higher than the prediction, but considering the error bars 
and the fact that the experimental ratio includes all kaons down to m± ~ 1 GeV, 
where the resonance contribution in particular to the K~ begins to become important, 
the agreement is good. 

The kaon to (anti-)lambda ratios, which are the only two ratios combining mesons 
with baryons, are quite badly off the predicted values; only for the lowest assumed 
freeze-out temperature the prediction is anywhere close to the data. This raises the 
question how strangeness neutrality is realized in the experimental data. Clearly more 
work is required to clarify these issues. 

There is also a recent experimental value for A/p for the similar S+Au collision 
system from the NA35 coUaboration-^^^ at CERN. Unfortunately, it was measured in 
a different kinematic domain, namely forward of central rapidity (3 < y < 5) and 
integrated over all p± < 2 GeV. Therefore the comparison with WA85 is somewhat 
problematic. The experimental value from NA35 is A/p = 0.8 it 0.3. The ratios 
predicted from the thermal parameters in Table 1 for the same p± range are 0.59ib0.06 
(A), 0.60 lb 0.06 (B) and 0.57 it 0.06 (C), respectively, for the three scenarios. Again 
there is agreement, but due to the rather large error bars and the different y- window 
of the data the value of this statement is limited. 

From this comparison, taking into account the partially still very large experi- 
mental uncertainties, we conclude that the concept of relative chemical equilibrium 
works reasonably well within the baryonic and mesonic sectors separately, but fails 
badly when mesonic and baryonic yields are compared with each other. We will en- 
counter a similar problem with pions relative to the baryons when we analyse entropy 
production in Section 7. This indicates a breaking of chemical equilibrium between 
mesons and baryons. 
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6.2 Thermal and Chemical Parameters from the NA35 Experiment 

For the strange particle production data of the NA35 collaboration a similar analysis 
was done in Ref. 20, but now for the lighter, symmetric system S+S. Although this 
experiment has not measured any multistrange baryons or antibaryons, it has accu- 
mulated a complete set of m_\_ and y-spectra for all the hadrons which they could 
identify. Due to the symmetry of the system, data taken at rapidities y < 3 can 
be simply reflected around the nucleon-nucleon center of mass at = 3 to obtain 
a complete rapidity spectrum. Therefore this experiment can be used to perform a 
global chemical analysis using the phase space integrated particle yields (47r data), or 
to study the rapidity dependence of the chemical parameters. 

The advantage of using full phase space data lies in the insensitivity to the flow 
dynamics at freeze-out. In the WA85 analysis longitudinal and transverse flow com- 
ponents had to be taken into account by cutting the theoretical flow spectra to the 
experimental kinematic acceptance range. For the A.tt data this is not necessary. How- 
ever, to achieve this simplification, one must assume that the system is characterized 
by a single set of thermal freeze-out parameters. While the observed rapidity in- 
dependence of the slope of the m_\_ spectra suggests a roughly rapidity independent 
freeze-out temperature and transverse flow velocity, the assumption of rapidity in- 
dependent chemical potentials cannot be justified. A A.tt analysis can thus only give 
information on the rapidity averaged chemical potentials. 

In Ref. 20 a chemical analysis of the A.tt data was performed and compared with an 
analysis of the particle ratios from the midrapidity bin only. Except for the extracted 
values for the baryon chemical potential (whose A.tt average came out larger than the 
central rapidity value, in agreement with expectations based on the central rapidity 
dip of the proton distribution seen in this collision system), the results and conclusion 
were the same. More recently, Slotta et al. performed a fit to the complete rapidity 
dependence of all particle yields and ratios. We will comment on that work at the 
end of this subsection and show that it mostly confirms the conclusions of Ref. 20. 

Within the global ansatz for the 47r-integrated data one can also easily compare 
the total hadron production from S+S and from elementary nucleon-nucleon (N+N) 
collisions at the same collision energy per nucleon. Gazdzicki et al.^°^ have extracted 
total multiplicities for an isospin symmetric nucleon-nucleon collision at 200 A GeV by 
isospin averaging available p+p and p+n (p+d) results, using a reasonable extrapola- 
tion to the n+n case. A comparison with the isospin symmetric S+S collision system 
illustrates very clearly the genuine nuclear (medium) effects on particle production. 
We will perform such a comparison within the thermal model and parametrize these 
effects via the thermal and chemical parameters. In spite of all the caveats about 
applying a thermal model to single N+N collisions, the result provides a very nice 
intuitive understanding of the basic difference between nucleon-nucleon and nucleus- 
nucleus collisions. 

The relevant data"^^"^'"^^^'"^^^""^^^ are summarized in Table 3. From the quoted 
multiplicities we build three ratios, avoiding meson to baryon ratios because of the 
indications for a meson/baryon non-equilibrium from the WA85 analysis above. 
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Table 3: List of measured particle multiplicities and corresponding ratios. The S+S data are from 
Refs. 102, 157-165, the N+N data from Ref. 101. The values listed under "S+S midrapidity" are 
the hadron multiplicities in the rapidity interval 2 < y < 3. 



Multiplicity 
& Ratio 


S+S 


s±s 

midrap. 


N±N 


K+ 
K- 

A 

A 

V — V 
N- 


12.5 ±0.4 
6.9 ±0.4 
8.2 ±0.9 
1.5 ±0.4 
20 ± 3 
98 ±5 


3.2 ±0.5 
2.2 ±0.5 
2.05 ±0.2 
0.57 ±0.2 
3.2 ± 1 
26 ± 1 


0.24 ±0.02 
0.17 ±0.02 
0.096 ±0.015 
0.013 ±0.005 
0.90 ± 0.09 
3.22 ±0.06 


o _ A 

Rk- 
o _ A 

-Ka.p — 

p-p 


0.18 ±0.05 

1.8 ±0.1 
0.41 ±0.08 


0.28 ±0.1 
1.45 ± 0.4 
0.64 ±0.2 


0.135 ±0.055 
1.4 ±0.2 
0.11 ±0.02 


Da 


0.09 ±0.02 


0.06 ±0.02 


0.12 ±0.01 



First estimates for NA35 

We begin again by using the approximation (75) to obtain a first estimate of the 
thermal parameters. From the two relations 

i2A = A-nr% Rk = KK" (92) 

one easily gets 

Aq = Rf^RT = 1-47 ± 0.02 , A3 = Rl^'^R'^'^ = 1.09 ± 0.02 . (93) 

With the available data 75 cannot be estimated by simple algebra. For this we would 
need the ratio -Ra.p = A/p in order to combine it with -Ra.p such that the rest masses 
and the temperature cancel. This ratio was recently published-^^^, but only in a 
restricted kinematic region. Therefore we must determine 75 numerically from a 
complete resonance gas analysis. 
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Resonance gas analysis for NA35 

The treatment of particle production in an equilibrated hadron gas is always sensitive 
to its mass spectrum {i : gi^rrii}. It should contain all independent hadronic degrees 
of freedom at freeze-out. However, resonance data become sparse above 1.8 GeV, 
and thus we don't have a complete set of hadronic states. In contrast to the WA85 
analysis, where we always used all well known^^ resonance states up to a mass of 2 
GeV, we will here also explore the sensitivity of the chemical analysis to this limit 
of our knowledge by restricting the experimental mass spectrum in alternative ways, 
labeled by A, B and C(mcut): 

A: Only the measured particles tt, if, p, (n). A, plus the S are taken into account; 
the S is added because a large fraction of the measured A stems from S° — )• A+7. 

B: We include the pseudoscalar and pseudovector meson nonet and the baryon 
spin- 1/2 octet and spin-3/2 decuplet. 

C: We include the full particle spectrum of all known hadronic resonances with 

We also treat the temperature differently than before: Whereas in the WA85 
analysis it was essentially taken as a free parameter, we fix it here from the condition 
of strangeness neutrality, since certainly within the full phase space strangeness must 
be conserved. 

To obtain the parameters T, Aq, Ag and 75, we now must solve a similar set 
of equations as in (89)-(91) plus the strangeness neutrality condition (56). The 
contributions from strong and electromagnetic decays after freeze-out are taken into 
account. Decay products from weak decays can be neglected because they are strongly 
discriminated against in the experiment by the requirement that the reconstructed 
tracks point back to the target. The results are given in Table 4. Although for 
the measured particle species the resonance decay contributions are appreciable, the 
simple first estimates (93) are again not strongly affected, because of a large amount 
of cancellation of these contributions from particle/antiparticle ratios. 7s is mostly 
sensitive to the experimental A/(p — p) ratio for which no such cancellation occurs. 
Accordingly we see a strong sensitivity of 7s to the chosen resonance spectrum. This 
is also graphically presented in Figure 12. The errors on the freeze-out parameters 
were calculated by error propagation from the measured multiplicities. 

We now shortly discuss these results with respect to their physical and systematic 
reliability. The physical validity of our procedure can be tested by comparing these 
full phase space results with a similar analysis of midrapidity ratios. Of course, the 
midrapidity results are correlated to the full phase space results since the same data 
set is used; an estimate of the systematic uncertainties of our procedure can, however, 
still be obtained from this comparison. Since the midrapidity ratios have much larger 
statistical errors, error propagation becomes a problem, and we thus refrain from 
giving error estimates for the midrapidity results in Table 4. 
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Table 4: Generalized thermal model parameters resulting from our analysis of the particle ratios 
given in Table 3, for the various tested scenarios as described in the text. 











7s 


T 
1 

[MeV] 




A 


1.47 ±0.07 


1.09 ±0.05 


0.58 ±0.13 


191 ± 13 




B 


1.48 ±0.07 


1.08 ±0.06 


0.67 ±0.15 


187 ± 16 


s+s 


C(1.4) 


1.51 ±0.07 


1.04 ±0.07 


0.64 ±0.14 


217 ±30 




C(1.6) 


1.52 ±0.08 


1.02 ±0.08 


0.84 ±0.18 


219 ±35 




C(1.8) 


1.54 ±0.06 


1.00 ±0.09 


0.95 ±0.20 


205 ± 39 




C(2.0) 


1.57 ±0.06 


0.99 ±0.09 


1.00 ±0.21 


197 ± 29 




A 


1.32 


1.09 


0.76 


185 




B 


1.32 


1.08 


0.88 


176 


s+s 


C(1.4) 


1.33 


1.06 


0.83 


198 


midrap. 


C(1.6) 


1.34 


1.05 


1.05 


193 




C(1.8) 


1.35 


1.04 


1.16 


183 




C(2.0) 


1.36 


1.04 


1.19 


178 




A 


1.48 ± 0.12 


1.25 ± 0.12 


0.17 ±0.04 


154 ± 22 




B 


1.48 ± 0.12 


1.25 ± 0.11 


0.19 ±0.04 


152 ± 21 


N+N 


C(1.4) 


1.51 ± 0.12 


1.18 ±0.16 


0.18 ±0.04 


174 ± 43 




C(1.6) 


1.52 ± 0.12 


1.18 ±0.15 


0.21 ±0.05 


172 ± 44 




C(1.8) 


1.51 ± 0.12 


1.18 ±0.16 


0.21 ±0.05 


163 ± 35 




C(2.0) 


1.51 ± 0.12 


1.18 ±0.16 


0.22 ±0.05 


161 ±31 



The observed differences to the full phase space results are consistent with simple 
expectations: from the shape of the proton rapidity spectra-^^^ it is known that the 
rather small sulphur nuclei cannot fully stop each other, resulting in a baryon number 
deficiency at midrapidity. On the other hand, excess strangeness production was 
found to be maximal near midrapidity^^^"^^^. Correspondingly we find a reduced 
value of Aq and an increased value of 7s at midrapidity. Values of 7s above unity 
should not be taken too seriously: the use of the strangeness neutrality condition 
in a restricted rapidity window is not absolutely safe, and the large experimental 
error (30%) on the experimental A/(p — p) ratio also directly propagates into the 
determination of 7s. 

Altogether the differences between the midrapidity and the 47r values are small, 
and for As and T the differences are even below 10%. This was recently confirmed in 
an analysis by Slotta^^ of the full rapidity spectra from this experiment. Assuming a 
boost-invariant longitudinal expansion velocity profile with a maximum flow rapidity 
ff of ±1.75 relative to the center of mass (this value yields an excellent fit to the shape 
of the measured pion rapidity spectrum), a rapidity independent freeze- out temper- 
ature T and strange quark chemical potential /is = 0, he was able to reproduce the 
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Figure 12: The dependence of the extracted thermal parameters Aq, As, T and 7s on the mass 
spectrum of the hadron system at freeze-out for central S-S collisions at 200 A GeV (47r data). (For 
details see text.) Also one set of values for N+N collisions is plotted. 



Pt and rapidity spectra of the measured proton and all strange hadron spectra from 
the NA35 S+S experiment simultaneously, in shape as well as normalization. The 
calculation contains three free parameters: the value of the baryon chemical potential 
at midrapidity /ig, one parameter A describing its dependence on the longitudinal 
coordinate (resp. the space-time rapidity ?/), and a strangeness suppression factor 
7s. /ig and A were fixed from the relative normalization and shape of the A and A 
rapidity spectrum using the parametrization fiB{v) = /^b + where rj is measured 
from midrapidity. 7s was determined by adjusting the relative normalization of the 
proton and A spectra. Assuming a freeze-out temperature of T = 200 MeV, Slotta 
found n% = 177 ± 33 MeV, A = 34.5 ± 1.5 MeV, and 7s = 0.75 (see Figure 13). 
Averaged over rapidity, the baryon chemical potential comes out as 240 it 40 MeV, in 
good agreement with the 47r analysis which gives 264 it 72 MeV (see Table 4). Over- 
all strangeness neutrality is satisfied at the level of a few percent; perfect neutrality 
could be reached by a slight adjustment of the temperature. Problems arise with the 
normalization of the pion spectrum which is underpredicted by about a factor 2 - 
this problem is generic and will be discussed extensively later in this review. 

In Figure 14 an attempt is shown to fit the same rapidity spectra with a lower 
freeze-out temperature of T = 150 MeV, compensating the for the slope of the exper- 
imental m± spectra (which for the S+S data corresponds to an apparent temperature 
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yD U \ d ^ ^ O y O 

Figure 13: Rapidity distributions of A, A, K+, K", K° and p - p from the NA35 S+S 
experiment^°^'^^^'^^^. The theoretical curves assume T = 200 MeV, rj = 1.75, and /is = ev- 
erywhere. The other parameters are described in the text. The theoretical spectra are cut to the 
experimental acceptance windows, indicated in the figures by the index "e" on the vertical axis. 
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of about 200 MeV) by a transverse flow velocity f3{ = 0.31. This is in the spirit of sce- 
nario B of the WA85 analysis. The rapidity dependent baryon chemical potential is 
now given by /ig ~ 130 MeV and yl ~ 24 MeV. In this case the proton spectra are au- 
tomatically correctly normalized with 7s=l, but the strangeness neutrality and kaon 
normalization are badly off - very similar to the observations made in the WA85 S+W 
case. It is striking to see that one additional parameter, namely a meson/baryon sup- 
pression factor = 0.37, repairs both^^ deficiencies very economically (Figs. 14d,e). 
The problem of underpredicting the pions^^ by a factor 2-3 remains, however. 

Let us now return to the sensitivity of the chemical analysis to the hadron mass 
spectrum. This is illustrated in Figure 12 for the analysis of the 47r data of NA35. 
From left to right the assumed spectrum of hadron resonances gets more and more 
rich. While T, Aq and Ag appear to be quite insensitive to this, 7s shows appreciable 
variations with the hadron mass spectrum. However, as more and more known reso- 
nances are included in the spectrum, all values seem to converge to asymptotic limits 
given by As ~ 1, Aq ~ 1.6, 7s ~ 1 and T ~ 200 MeV. 

We conclude that chemical freeze-out in central S+S collisions at 200 A GeV is 
characterized by: 

• a vanishing strange quark chemical potential (As ~ 1 or /is ~ 0); 

• a largely saturated strange quark phase space (7s ~ 0.75 — 1). 

The chemical freeze-out temperature is around T = 200 MeV if relative chemical 
equilibrium and strangeness neutrality of the hadron gas are assumed. However, we 
argued before that this temperature is uncomfortably high. Furthermore, the hadron 
gas picture cannot correctly reproduce the pion multiplicity. Reducing the freeze-out 
temperature to more "reasonable" values by postulating transverse flow leads to a 
violation of the strangeness neutrality. This can be repaired by assuming a global 
meson/baryon suppression (i.e. by giving up relative chemical equilibrium between 
mesons and baryons), but this does not resolve the pion deficiency problem. All of 
these features are in qualitative and quantitative agreement between the S+W and 
S+S collision systems. 

An interesting additional result is the difference between the S+S and the N+N 
data, shown in the right sector of Figure 4 for the C(2.0) scenario. Although one 
must be careful when applying the thermal description to N+N collisions, we think 
that the extracted results still show clearly the important qualitative difference be- 
tween nucleon-nucleon and nucleus-nucleus collisions: While Aq takes similar values, 
intermediate between the S+S 47r and midrapidity results, the strange quark fugac- 
ity tends to deviate from unity. The temperature is somewhat lower than in S+S 
case (it is close to the canonical Hagedorn value of 160 MeV); this is also seen in 
the slope of the m± spectra. A clear qualitative difference is observed, however, for 
the strangeness saturation factor: it is 4 times larger in central S+S collisions than 
in N+N interactions. The strangeness enhancement is therefore clearly seen in the 
thermal description by an increased 7s factor. 
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Figure 14: Similar to Figure 13, but 
now assuming T = 150 MeV together 
with a transverse flow velocity /3f = 
0.3. In this case the protons are cor- 
rectly reproduced with 7s = 1, but 
the correct kaon normalization and over- 
all strangeness neutrality require a me- 
son/baryon suppression factor jy^ = 
0.37. 
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Table 5: Predicted total particle multiplicities in 200 A GeV S+S collisions for scenario C(2.0), 
with the freeze-out volume fixed by the total K'^ multiplicity. In the lower part of the table we 
give the abundances of hadronic resonances before their strong decays, but including the feed-down 
contributions from higher resonances. 



Total 


S+S 


S+S 


N+N 


multiplicities 


Att 


midrap. 


Att 


K+ 


12.5 


3.2 


0.24 


N- 


65 


15 


3.9 


7r~ 


54 


12 


3.6 


K- 


6.9 


2.2 


0.17 


A 


7.5 


1.4 


0.097 


A 


1.35 


0.39 


0.013 


P-P 


18 


2.2 


0.89 


P 


1.4 


0.43 


0.081 




0.73 


0.19 


2.4 X 10"^ 




0.31 


0.089 


5.5 X 10"^ 


n 


0.11 


0.037 


7.5 X 10-5 


n 


0.12 


0.029 


2.8 X 10-5 


V 


4.9 


1.29 


0.34 


/ 


6.7 


1.48 


0.48 




5.5 


1.32 


0.45 


</> 


1.3 


0.43 


5 X 10-^ 


A 


17 


2.2 


0.84 


A 


1.17 


0.35 


0.071 



Predictions resulting from NA35 

Similar to the WA85 case we can use the extracted thermal parameters to predict 
further particle ratios^^. Table 5 lists the predicted total multiplicities from which 
any ratio can be constructed. The (point particle) volume parameter (46) was 
fixed by the total yield and then used to normalize all other multiplicities. 

Comparing the predictions of Table 5 with the measured values of Table 3 the 
following points emerge: 

• The K~ yield is exact, because the yield and the if"*" /if- -ratio were fitted. 

• For the baryons A, A and p — p data and prediction agree within error bars; 
this indicates that the K/ A ratio is correctly reproduced - all other ratios were 
already used in the thermal analysis. 

• Table 5 gives a A/p ratio of 1.0 (47r) and 0.9 (midrapidity). The recently 
published-^^^ central rapidity value from the NA35 collaboration of A/p = 1.5 ± 
0.4 is somewhat on the high side, but also affected by large statistical errors. 
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Table 6: Particle ratios and acceptance windows covered by the NA36 S+Pb experiment at CERN. 
The ratios in the lower part of the table were obtained by extrapolating the data to a common 
mx-window using a thermal ansatz (59) with the measured mx-slopes, but without correction for 
the different y acceptances. All values are taken from Ref. 173. 



Ratio 


value 


y range 


px -range 


Ra = A/A 
Re = H+/S- 
Rea = S-/A 
%A = S+/A 


0.117 ±0.011 
0.276 ±0.108 
0.066 ±0.013 
0.127 ± 0.022 


2.0 < y < 2.5 
2.0 < y < 2.5 
1.5 < y < 2.5 
2.0 < y < 3.0 


0.6 < px < 1.6 GeV 
0.8 < px < 1.8 GeV 
0.8 < px < 1.8 GeV 
0.6 < px < 1.6 GeV 


Ratio 


value 


y range 


mx -range 


Ra = A/A 
Re = H+/S- 
Rea = S-/A 
Re,a = S+/A 


0.093 ±0.017 
0.276 ±0.108 
0.122 ± 0.024 
0.281 ±0.050 


2.0 < y < 2.5 
2.0 < y < 2.5 
1.5 < y < 2.5 
2.0 < y < 3.0 


1.5 < mx < 2.2 GeV 
1.5 < mx < 2.2 GeV 
1.5 < mx < 2.2 GeV 
1.5 < mx < 2.0 GeV 



• The negative particle multiplicity N is strongly underestimated. This will be 
discussed in Section 7 in the context of entropy production. 

• The resonance yields listed in Table 5 show the importance of resonance con- 
tributions. For example the four A isobars contribute about 50% to the final 
proton yield, and similarly for the antiprotons. 

6.3 Thermal and Chemical Parameters from the NA36 Experiment 

The NA36 collaboration has published-^^^'-^^^ strange baryon ratios from S±Pb col- 
lisions at 200 A GeV. They study the same set of ratios as the WA85 collaboration 
(Eqs. (79)-(82)), and since the target nuclei are not too different the results should in 
principle be comparable. Unfortunately, the NA36 ratios (listed in Table 6) are not 
all given in the same kinematic window, which makes the analysis more difficult. The 
rapidity spectra from S±S collisions in Figure 13 show that in general the particle 
ratios have a strong rapidity dependence, and since for the asymmetric S±W and 
S±Pb collisions the shape of the rapidity spectra is not known, a straightforward 
comparison between the two experiments is not possible. 

Although the ratios may be safe, we should still mention that the absolute nor- 
malization of the NA36 strange particle yields has been controversial (see NA35 and 
NA36 contributions in Ref. 12): for so far unknown reasons the A rapidity distri- 
butions from the NA35 and NA36 experiments disagree by about a factor 2 in nor- 
malization and also show a different shape. In view of these unsettled questions all 
theoretical interpretations of these data sets should be taken with a grain of salt. 
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First estimates for NA36 

Similar estimates as in Section 6.1 can be extracted from these data only if the rapidity 
dependence of the ratios is neglected. As already pointed out, this is very dangerous, 
and indeed a comparison of the results given in Refs. 172 and 173 shows that at least 
R\ is very sensitive to the selected rapidity interval: 

1.5<y<3.0: i^A = 0.207 ± 0.012 (Ref. 172) ; 

2.0<y<2.5: i^A = 0.117 ± 0.011 (Ref. 173) . (94) 

The p± range for both ratios is the same (0.6 < p_L < 1.6 GeV). 

Since the particle-antiparticle ratios in Table 6 were obtained in the same rapidity 
interval, one can at least apply Eqs. (85) and (86) to extract the fugacities-^^^ : 

Aq = 1.78 ±0.09; A3 = 1.03 ± 0.06 . (95) 

The estimate for 75 from Eq. (87) is less reliable, due to the different rapidity intervals 
of the required ratios. Indeed, the identity (-Ra/-Rh)(-Rh a/-Rh,a) = 1 is violated 
by the ratios given in Table 6 by 25%. Neglecting this problem one obtains 75 = 
0.19 lb 0.02, and after multiplication by a factor 2 for the S° contribution one gets^^^ 
7s = 0.38 lb 0.04. This is below the naive estimate (88) from the WA85 data, but we 
already know that it will be corrected upward once resonance decays are included. 



Resonance gas analysis for NA36 

The NA36 coUaboration-^^^ have analyzed their data following the work of Cleymans-^^ 
et al. This incorporates the resonance contributions. In this description strangeness 
neutrality is taken as a constraint, and full chemical equilibration, 7s = 1, is assumed. 
The result is shown in Figure 15. All ratios are compatible with a fully equilibrated 
hadron resonance gas at T = 172 ib 16 MeV and fiB = 290 ib 50 MeV. This result can 
be summarized by the parameters 

Aq = 1.76 ±0.15; As = 1.07 ± 0.07 ; (96) 
r = 172 ±16 MeV; 7s = 1.0; (97) 

After the first naive estimate for 7s in Section 6.3 came out below the corresponding 
one for the WA85 results, it is surprising to see that this resonance gas analysis 
yields compatibility of the NA36 data with 7s = 1, while a similar analysis of the 
WA85 data-^^ finds strong incompatibility with 7s = 1 and rather confirms the value 
7s ~ 0.7 extracted in Section 6.1. On the other hand, the resulting values for Aq and 
T are also different. These differences are obviously a consequence of the different 
experimental particle ratios, in particular of R\ which differs from the WA85 value by 
more than a factor 2. As long as quantitative experimental information on the crucial 
rapidity dependence of these ratios is not available, more detailed statements about 
the consistency of the two experiments and their respective theoretical interpretations 
cannot be made. It is, however, interesting to observe that the strange quark chemical 
potential seems to vanish in all these analyses of heavy ion data at CERN energies, 
independent of the selected rapidity window, and 7s is always quite large. 



54 



Figure 15: The T — /iB-plane for a 
strangeness neutral resonance gas with 
7s = 1. The values compatible with the 
measured particle ratios within one stan- 
dard deviation (after resonance decays) are 
shown as hatched areas. The common 
overlap region defines the boundary of al- 
lowed values for T and ji^ ■ The Figure was 
taken from Ref. 173. 



6.4 Thermal and Chemical Parameters from the AGS Experiments 

The heavy ion experiments at the AGS in Brookhaven are performed at a lower beam 
momentum of about 10 - 15 A GeV. Comparing these experiments with the ones at 
CERN energies of 160 - 200 A GeV allows to study the beam energy dependence 
of particle production. We are interested in the differences in strange hadron pro- 
duction and in the resulting changes in the thermal fireball parameters. An analysis 
along these lines of the available data at the AGS from Si+Au (Pb) collisions was 
recently done in Ref. 26. Let us shortly discuss some small differences relative to the 
investigation of the CERN data presented above: 

• The hadron gas used in Ref. 26 incorporates baryons up to masses of 2 GeV, 
but mesons only up to 1.5 GeV. In view of the lower temperature at the AGS 
(see below) this seems justified. 

• In Ref. 26 a finite size correction-^^^ to the density of states in the partition 
function was applied. This correction is of the order ~ Ath/-R where Ath = 

■j2'K/mT is the thermal wavelength of a particle with mass m and R is the linear 
dimension of the system. This correction is important for low temperatures 
and for light particles and mostly affects the ratios between a heavy and a light 
particle (e.g. 7r/(p + n)). 

• For extensive observables the volume correction was taken into account via the 
excluded volume description^^. This drops out from particle ratios. 

• The calculation was done assuming absolute chemical equilibrium (75 = 1) from 
the beginning, i.e. no strangeness suppression was allowed for. 
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Table 7: Particle ratios calculated in a thermal model for two different temperatures, a baryon 
chemical potential of /iB = 0.54 GeV and a strange quark chemical potential jig fixed by the condition 
of strangeness neutrality, in comparison with experimental data (statistical errors in parentheses) 
from central Si + Au (Pb) collisions at 14.6 A GeV. The table was taken from Ref. 26. 



Particles 


Thermal Model 
1 =izU iVieV 1 =i4U iVieV 


Data 

exp. ratio rapidity Ref. 


7r/fp+n) 
d/(p+n) 
p/p 


1.29 1.34 

4.3 • 10-2 5.8 • 10-2 

1.47 • 10-4 5.8 • 10-4 


1.05(0.05) 0.6 - 2.8 139 
3.0(0.3) • 10-2 0.4 - 1.6 139 
4.5(0.5) • 10-4 0.8 - 2.2 149 


K+/7r+ 
K'/tt- 
K°/7r+ 
K+/K- 


0.23 0.27 
5.0 • 10-2 6.2 • 10-2 
0.14 0.16 
4.6 4.3 


0.19(0.02) 0.6 - 2.2 139 
3.5(0.5) • 10-2 0.6 - 2.3 139 
9.7(1.5) • 10-2 2.0 - 3.5 144 
4.4(0.4) 0.7 - 2.3 139 


A/(p+n) 
A/A 


9.5 • 10-2 0.11 

8.8 • 10-4 3.7 • 10-3 


8.0(1.6) • 10-2 1.4 - 2.9 144 
2.0(0.8) • 10-3 1.2 - 1.7 149 


</>/(K++K-) 


2.4 • 10-2 3.6 • 10-2 


1.34(0.36) • 10-2 1.2 - 2.0 149 


E-/A 


6.4 • 10-2 7.2 • 10-2 


0.12(0.02) 1.4 - 2.9 146 


d/p 


1.1 • 10-^ 4.7 • 10-^ 


1.0(0.5)- 10-5 2.0 108 



Together with the requirement of strangeness neutrality, this last assumption re- 
duces the number of input parameters to two: Tf and fiB (/^q)- The temperature 
was fixed by an analysis of the measured A(1232) contribution to the proton yield. 
An investigation-^^'* showed that the A/p-ratio is compatible with a temperature in 
the range of T = 120-140 MeV. This very useful information is, unfortunately, not 
available from the CERN experiments. If measured well, it allows an accurate extrac- 
tion of the chemical freeze-out temperature independent of the baryon and strange 
chemical potentials. 

The baryon chemical potential is fixed by a compromise between the measured 
7r/(p + n)-ratio and considerations regarding the particle density at freeze-out^^. The 
latter can be fixed from the size of the freeze-out volume as measured by HBT 
interferometry-^^^ and the measured baryon multiplicity. The resulting value was 
given as^^ /xb = 540 MeV. 

The results of this analysis are shown in Table 7. In general the agreement be- 
tween the thermal model ratios and the experimental ones is very good. The an- 
tibaryon/baryon ratios are very sensitive to the temperature for fixed fiB, in the given 
temperature interval they vary by a factor of 4 and cover a range which includes the 
experimental value. 

In view of the different rapidity windows for the various measurements and the 
generic rapidity dependence of particle ratios such an agreement is not trivial. How- 
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ever, most of the data cover a rather wide rapidity range around the central rapidity"^ 
of 1.4, and stopping is more efficient and longitudinal flow less important at AGS 
energies-^^^'-^^^. To a large extent the analysis of Ref. 26 thus resembles in character 
the global analysis of the 47r NA35 data at CERN in Section 6.2. 

Taking a closer look at the strange particle ratios one observes an interesting sys- 
tematic behavior. While within the chosen temperature interval the predictions for 
the two particle/antiparticle ratios K'^ j K~ and A/ A overlap nicely with the exper- 
imental values inside the error bars, this is not the case for some of the other ratios 
containing strange quarks. Let us therefore look at those and divide the experimental 
mean values by the two model values. Neglecting the S/A-ratio (see below) we obtain 

ratio T = 120 MeV T = 140 MeV 



K+/Tr+: 0.82 0.70 

exp. value ^ K'/tt-: 0.70 0.56 

model value K°/Tr+: 0.70 0.61 

A/{p + n): 0.84 0.72 

^/{K+ + K-): 0.56 0.37 



The model ratios are systematically below the experimental values by a factor 0.6-0.8. 
Now all these ratios involve particles with one more strange quark or antiquark in the 
numerator than in the denominator. This indicates incomplete strangeness satura- 
tion and thus a violation of the model assumption of absolute chemical equilibrium. 
We can soften this assumption by assuming only relative chemical equilibrium and 
introducing as before the parameter 7s through Eq. (38). From the above table we 
extract a value of 7s = 0.7 it 0.1; thus the level of strangeness saturation reached at 
the AGS is similar to the CERN value. 

We left out the S/A-ratio for a particular reason: neither does it follow the above 
systematics, nor does it agree with the predictions of any of the available event 
generators''*^. So far there is no model which can explain such a large ratio con- 
sistently with the other observations. This strange result cries out for independent 
experimental confirmation. 

Braun-Munzinger et al.^^ have also tried to confirm the thermal parameters by 
comparing the measured transverse momentum spectra with the model predictions 
from a thermalized source with the same temperature undergoing longitudinal and 
transverse hydrodynamical fiow^^'^^. Since the model was developed for symmetric 
systems and the new momentum spectra for Au+Au collisions-^^^'-^'*^ were not yet 
publicly available, this comparison was done in Ref. 26 for the lighter Si+Al system^. 
The mj^-spectra were described by Eq. (60) with a = 1. After adding resonance decay 
contributions excellent agreement^^ is obtained for all particle spectra from the pions 
to the deuterons. The presence of transverse fiow is very clearly seen in the slope of 

"^We define as central rapidity the value where the pion production peaks-^^^. 

similar comparison of the m± spectra with E802 data from Si+Au collisions is shown in 
Ref. 107, with very similar parameters. 
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the experimental spectra for m± — mo > 0.3 GeV (i.e. beyond the resonance decay 
region) which flattens considerably with increasing mass of the hadrons^^; in the data 
presented in Ref. 26 the deuteron and pion slopes differ by about 50%. This is well 
reproduced by the assumed average flow velocity (/3f) = 0.39 (0.33) for Tf = 120 (140) 
MeV, respectively. 

We conclude that a thermal model for Si+Au(Pb,Al) collisions at AGS energies 
works well and leads to the following thermal parameters: 

T = 130 ± 10 MeV /ib ~ 540 MeV . . 

7s = 0.7 ±0.1 /is = 121 ± 14 MeV ^ ^ 

In terms of the fugacities used before this translates as 

Aq = 4.0±0.5 ; As = 1.55 ±0.1. (99) 

This result can be compared with the earlier, less elaborate thermal model in- 
terpretations of Refs. 23, 24. The authors of Ref. 23 consider only the two parti- 
cle/antiparticle ratios K'^ /K~ and A/ A from Table 7 as input. As discussed before, 
resonance decay contributions can be neglected for these ratios. Using the require- 
ment of strangeness neutrality the fugacities and also the freeze-out temperature can 
be extracted^'*: 

Aq = 3.7±0.3 ; As = 1.72 ± 0.19 ; Tf = 127 ± 8 MeV . (100) 

This agrees very nicely with Eqs. (98), (99); the small differences are probably due 
to the fact that, in contrast to Ref. 23, the A/A-ratio was not fitted exactly in Ref. 26. 



7 Specific Entropy from Heavy Ion Experiments 

The careful reader will have noticed the conspicuous absence of any discussion of 
the pion yields so far. Although the pion is the most abundantly produced hadron, 
especially so at CERN energies, it carries "no" conserved quantum numbers and is 
thus unsuitable for extracting information about the chemical potentials. In chemical 
equilibrium the pion multiplicity is a function of the temperature alone. However, of 
all the thermodynamic parameters used in Section 6, the chemical freeze-out temper- 
ature Tf is the most uncertain one. This uncertainty affects the prediction of total 
pion multiplicities in a very direct way. 

On the other hand the pions, being the most abundant of all produced hadrons, are 
the dominant carriers of entropy in the collision fireball. This was already discussed 
in Section 5.5 where we related the intensive quantity S/ B to the specific net charge 
Dq which is normalized to the total charged particle multiplicity and thus dominated 
by pion production. 

In Section 6 we discussed extensively whether at freeze-out the collision zone can 
be described as a hadron resonance gas in thermal and some kind of relative chemical 
equilibrium. We will now discuss and test an additional aspect of this picture, namely 
its specific entropy content. 
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A measurement oi Dq as a function of rapidity was performed for the S+Pb system 
at 200 A GeV by the EMU-05 collaboration^"^ The value at central rapidity is Dq = 
0.085 lb 0.01. This can be compared with values calculated from the hadron gas model 
with the thermal parameters extracted from the WA85 data. (Both collision systems 
feature nearly the same number of target participants.) No detailed comparison has 
been made for the NA36 parameter values, but since the thermal parameters extracted 
from this experiment are similar, the conclusions would be the same. 

The thermal model results for Dq and S/ B corresponding to the thermal param- 
eters extracted from the chemical analysis are given in the last row Table 1. Only 
for scenario B, which is plagued by a strong violation of strangeness neutrality, the 
model produces a value of Dq which is compatible with the experimental value of 
0.085 lb 0.01. In all other cases the calculated Dq is higher then the measured one, 
implying that the data contain a higher total charged multiplicity and larger specific 
entropy S/B than provided by the hadron gas model^°°. The discrepancy is of the 
order of a factor 2. 

The situation is similar for the NA35 experiments. The measured Dq is given in 
Table 3. The hadron gas values for the parameters from the scenario C(2.0) would 
be Dq = 0.15 for the 47r-data and Dq = 0.09 at central rapidity. In other words, the 
measured negative particle multiplicity [N~ = 98 ib 5 from Table 3) is higher than 
the predicted one [N~ = 65 from Table 5) - again the discrepancy is large, of the 
order of 50%. Similar conclusions have been reached earlier by Davidson et al.^^. 

An independent way to arrive at this statement has recently been presented by 
Gazdzicki^^^. He looked at the systematics of pion production in symmetric nuclear 
reactions at beam energies from 2 A GeV (BEVALAC) to 200 A GeV (CERN) and 
compared it with pion production in elementary nucleon-nucleon collisions in the 
same energy range-^-^-^. According to his analysis, one sees for BEVALAC and AGS 
energies a suppression of pion production in nuclear collisions relative to nucleon- 
nucleon interactions^^^: 

= - - -0-35 • (101) 

{Np) {Np)aa {Np)nn ^ ' 

Here (tt) is the mean pion multiplicity for symmetric nuclear collisions (A+A) or 
nucleon-nucleon reactions (N+N), respectively, and {Np) is the number of partici- 
pating nucleons. Note that in the studied energy range (101) turns out to be valid 
independent of the center of mass energy ^/iiviv- Proceeding to the higher CERN 
energies, he obtains from an analysis of the NA35 data on S+S coUisions-^-^^ a value 
of 

(tt) 

A-^ = +0.9 ±0.6. (102) 

{Np) 

Gazdzicki also derives an absolute value for the entropy in the form S' j {Np) = 
6.8 lb 0.6, where S' is the entropy in units of the entropy per pion-^-^^. Taking for 
the latter the thermal value of about 4 units per pion, we obtain S j B = 27 ib 2. 
This is again about 50% higher than the hadron gas value resulting from the thermal 
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model parameters for the S+S collisions given in Table 4, in agreement with our above 
conclusion. 

This feature is generic for all heavy ion experiments at CERN energies. The 
hadron gas model, assuming relative chemical equilibrium with only the absolute 
strangeness level out of equilibrium, is not consistent with all experimental observa- 
tions. As long as the model is forced to yield zero net strangeness, it predicts values 
of Dq well above the measured values, irrespective of the particular choice for the 
freeze-out temperature^^. In other words, the specific entropy S j B in the CERN data 
is large and cannot be reproduced by the hadron gas model. 

At the lower AGS energies Dq has not been investigated experimentally. But we 
can consider the absolute pion production directly. Table 7 shows that the thermal 
model slightly overestimates the measured pion/nucleon ratio, but not at a level which 
would indicate serious disagreement. From the extracted thermal model parameters 
a value of S j B = 13 it 1 was calculated in Ref. 24; this is significantly lower than 
the CERN values. We conclude that at AGS energies the entropy production is low 
and can be described by the thermal model with a hadron resonance gas equation of 
state, while at CERN energies entropy production is dramatically higher and cannot 
be reproduced by the hadron gas model. 

8 Discussion of the Results 

We have shown the results of an interpretation of experimental particle abundances 
in a thermal hadron gas model, with special emphasis on strange particles. For the 
lower AGS energies the model yields an excellent description of all hadron spectra 
and particle abundances from Si+Al and Si+Au collisions. It will be interesting to 
see whether this can be confirmed for Au+Au collisions where some of the necessary 
data have already been taken; this largest presently available collision system at the 
AGS is expected to show even clearer signs of thermal and chemical equilibration, 
superimposed by collective expansion. 

At the higher CERN energies, on the other hand, some features of the data can 
not be reproduced by the model. What exactly goes wrong is not yet quite clear; the 
interpretation of how the model fails depends to some extent on how it is applied to 
the data. The crucial parameter is the freeze-out temperature whose determination is 
a bottleneck in the analysis. Whereas the baryon and strange fugacities can be very 
precisely given by a few particle/antiparticle ratios (85), (86), the temperature is much 
harder to determine. One of the reasons is that one must distinguish between the 
thermal and the c/iemica/ freeze-out temperature. The ^/lerma/ freeze-out temperature 
can be extracted from the slope of the mj^-spectra, although the separation of thermal 
and collective motion is not quite model independent. Dynamical arguments^^ suggest 
thermal freeze-out 150 MeV. The chemical freeze-out temperature is 

reflected by the particle ratios; ratios between heavy (baryons) and a light hadrons 
(mesons) show the largest sensitivity to the freeze-out temperature. 

This sensitivity can, however, be exploited only if prior to freeze-out chemical 
equilibrium was established to begin with. Whether this is true or not is an important 
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dynamical question whose answer would provide crucial information about various 
time scales in nuclear collisions. We have argued that chemical equilibration happens 
in stages, and that the equilibration of the overall strangeness level and, in a hadronic 
environment, of the baryon+antibaryon<-)-meson processes should occur most slowly. 
If the latter processes have not equilibrated, the abovementioned meson/baryon ratios 
are useless for a determination of the freeze-out temperature. 

We have therefore based our chemical analysis on particle ratios taken entirely 
from within the baryonic or from within the mesonic sector. It turns out that the 
baryonic ratios by themselves do not fix the freeze-out temperature very well; if there 
is a preference at all, it is that the (not very accurately known) triple/double strange 
ratios prefer a high freeze-out temperature (of the order of the mj^-slope of 230 
MeV). 

Another (in principle) important constraint is the condition of overall strangeness 
neutrality for the decaying collision region. Unfortunately it can also be exploited only 
if the system is chemically equilibrated or if chemical equilibrium is broken in par- 
ticular, well-controlled ways which we have labeled by relative chemical equilibrium. 
Assuming relative chemical equilibrium with respect to all processes except strange 
pair production, the strange and non-strange baryon/antibaryon and meson/meson 
ratios (except for those involving pions) are compatible with strangeness neutrality in 
a hadron resonance gas if the freeze-out temperature is around 190 - 200 MeV. Lower 
temperatures are inconsistent with strangeness neutrality in a hadron gas unless the 
meson/baryon equilibrium is broken, but even by allowing for such a possibility the 
disagreement with the measured pion yields cannot be resolved. 

Summarizing these comments we conclude that the CERN data do not allow for a 
consistent and comprehensive interpretation within the hadron gas model, even after 
allowing for some simple patterns for breaking absolute chemical equilibrium. Still, 
the chemical parameters extracted from the analysis presented in the previous chap- 
ters show a number of strikingly consistent features which are worth being exposed in 
a systematic way. In spite of all the caveats above, these features provide an intuitive 
and internally consistent picture of some of the main characteristics of high energy 
nuclear collisions. 

A priori nuclear collisions are characterized by the two controllable parameters 
^/s, i.e. the center of mass energy per (participating) nucleon p air, and .Bpartj the 
number of participating projectile and target baryons, which is related to the size of 
target and projectile nuclei and the impact parameter. Accordingly we will subdivide 
the observables into two categories, those which scale dominantly with ^/s and those 
which show mostly a scaling with -Bpart- We call these different scaling laws "energy 
scaling" and "size scaling", respectively. 

Figure 16 summarizes graphically the thermal parameters extracted in Sections 6 
and 7. The initial fireball temperature is related to the initial energy density and 
thus expected to exhibit mostly an energy scaling; a weaker size scaling arises from 
the different stopping power of light and heavy nuclei which transforms beam energy 
into internal excitation energy. But what is the expected scaling of the freeze-out 
temperature Tf? 
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Figure 16: Summary of the thermal parameters extracted in Section 6. In the last row of 
diagrams, the triangles correspond to the value of S/B in a hadron gas with the thermal parameters 
shown in the first 4 rows, while the crosses correspond to "measured" values of this quantity. The 
values of S/B in the first three columns are from Ref. 110; note that the value plotted by a cross in 
the Si-Au column is really from Si-Al collisions'^''. The experimental S/B value in the last column 
is from the Dq measurement^''^ using Eq. (78). 



If we assume that Tf is determined by geometry and occurs when the mean free 
path Xp is of the order of the transverse radius R± of the fireball, 

R^T,) ^ Xp{T,) ^ (103) 
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where n is the fireball density and (crv) the thermal cross section (23) at freeze-out, 
then we expect Tf to scale with the inverse of the system size. The reason is that 
both the particle density and the thermal cross sections rise with temperature (see 
e.g. Figure 4). Eq. (103) can depend on the beam energy only through the baryon 
stopping power via the relationship between the fireball density n and Tf; its energy 
scaling is thus expected to be very weak. 

This expectation is in stark contrast to the behavior seen in Figure 16 where Tf 
shows a clear scaling with both increasing beam energy and increasing size of the 
system. Note that we plotted in Figure 16 the chemical freeze-out temperature as 
determined by the particle ratios and the strangeness neutrality condition (and not 
the measured slope parameter of the mj^-spectra). As seen, Tf follows the scaling be- 
havior expected for the initial temperature To as discussed above. Somehow the final 
particle ratios seem to have a memory of the initial conditions of the fireball. Naively 
this would seem to contradict the thermal interpretation because thermalization by 
assumption implies memory loss via an increase in entropy. 

However, in our discussion of thermal vs. chemical freeze-out we noted the rather 
strong evidence in the transverse momentum spectra for the presence of collective 
transverse flow at freeze-out. In the presence of flow the geometric freeze-out cri- 
terium (103) should be replaced by a dynamical one because collective expansion leads 
to dynamical decoupling before the geometric freeze-out criterium is satisfied^^'^^^. 
Hydrodynamical simulations suggest that in nuclear collisions thermal freeze-out is 
entirely dominated by the collective dynamics^^'^^. It is more than likely that the same 
is true for chemical freeze-out because the only difference between the two freeze-out 
criteria results from the substitution of the total by the corresponding inelastic par- 
tial cross section. Since higher initial energy densities generate larger hydrodynamic 
flow and larger hydrodynamic flow leads to earlier dynamical freeze-out at a higher 
freeze-out temperatures^^'^^, the qualitative behavior of Tf seen in Figure 16 can be 
understood in this picture. In this sense the development of collective flow helps the 
system to "remember" the initial conditions of the fireball in spite of thermalization. 

The two fugacities Ag and Aq fall off going from left to right in Figure 16. At 
the AGS the large value of Aq reflects the high baryon density caused by the nearly 
complete baryon stopping at these energies. A part, but not all of the strong decrease 
of Aq at CERN energies can be understood by the larger freeze-out temperatures there; 
on top of this trivial effect it reflects, however, the onset of partial transparency at 
CERN energies. Compared to the strong jump between the lower (first column) 
and the higher beam energies (second to fourth column), Aq varies very little with 
the system size. It should be noted, however, that the values in the second and 
third column are from 47r particle ratios, whereas the last column is from the WA85 
midrapidity ratios. The analysis of Ref. 85 showed that Aq depends on rapidity and, 
at CERN energies, is lower at midrapidity than in the nuclear fragmentation regions. 
Comparing the Aq values at midrapidity from S+S (NA35), S+W (WA85) and S+Pb 
collisions (NA36) as they were given in Section 6, one can see a rising tendency with 
the size of the nuclear target. This reflects again the increasing baryon stopping 
power of larger nuclei which counteracts the onset of nuclear transparency at CERN 
energies. 
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Hence Aq is a measure for the efficiency of baryon number stopping of different 
size nuclei at different energies. At the relativistic heavy-ion collider RHIC {^/s=200 
GeV) one expects Aq 1 at central rapidity, independent of the size of the colliding 
nuclei. 

In first approximation Ag must follow the behavior of Aq because of the strangeness 
neutrality relation. However, it is not all trivial that for all nuclear collisions at SPS 
energies the value of Ag is approximately 1, corresponding to a vanishing strange 
quark chemical potential. (In N+N collisions this seems to be different, but due to 
the large error bar a safe statement cannot be made.) Indeed, this appears to be the 
most stable result of the chemical analysis, and it is independent of the actual value 
of the freeze-out temperature and how the strangeness neutrality is actually realized 
between the baryons and mesons. We already commented before on the intriguing 
agreement of this value with the one prevalent in a strangeness neutral quark-gluon 
plasma. To what extent this may be more than a diabolical coincidence will be 
discussed in Section 9.2. 

The strangeness suppression factor 7s shown in the fourth row of Figure 16 demon- 
strates clearly that the physics of nuclear collisions is different from N+N collisions. 
For nuclear collisions 7s is close to 1, while in the N+N case is much smaller, 7s ~ 0.2. 
This size scaling of 7s was worked out in more detail in Ref. 22. The increase of 7s 
with -Bpart is a clear sign that strangeness equilibration is caused by secondary in- 
teractions. The large degree of strangeness saturation in nuclear collisions which is 
refiected by 7s < 1 is a strong indication that unusually efficient and rapid mechanisms 
for strangeness production are at work. 

In the second but last row of Figure 16 we have plotted the average transverse 
collective flow velocity at the point of ^/lerma/ freeze-out, i.e. at the point of decou- 
pling of the momentum spectra. For this we need to know the experimental slope of 
the m± spectra and the thermal freeze-out temperature. At the AGS the transverse 
flow and the thermal freeze-out temperature can be determined separately from the 
momentum spectra, because fortunately also momentum spectra of very heavy parti- 
cles like the deuteron^^ are available; still, the determination is not yet very accurate, 
and lower thermal freeze-out temperatures down to T ~ 100 MeV, compensated by 
correspondingly larger transverse flow velocities, still give acceptable fits to the m± 
spectra. The identity of the thermal and chemical freeze-out temperatures at the 
AGS is thus not rigorously proven, although it seems to work quite nicely. At CERN 
energies such a clear separation of thermal and collective motion from the spectra is 
not (yet) possible-'^. The transverse fiow must be determined with the help of a theo- 
retical model for the dynamics of the fireball and involves a freeze-out criterium. Such 
calculations are only available for the S+S case at SPS energies^^. These calculations 
suggest thermal freeze-out at T = 150 — 160 MeV with (/3f) = 0.25 — 0.3; this value 
has been indicated in the S+S column. For the N+N case collective hydrodynamical 
expansion is not expected, so we entered the value (/3f) = in the second column. 
The value in the last column was obtained by assuming (for lack of better knowledge) 

•^A first deuteron spectrum was presented by the NA44 collaboration in Ref. 113, but these data 
have not yet been subjected to a full thermal+flow analysis. 
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that thermal freeze-out occurs also in the heavier systems at T = 150 MeV (see sce- 
nario B in Section 6.1), by averaging over the mj^-slopes quoted by WA85"^^^'"^^^ and 
NA36i^^ 

Generically, the experimental m± spectra flatten for larger collision systems, in- 
dicating earlier thermal decoupling at a higher temperature and/or larger collective 
transverse flow. Hydrodynamical calculations^^ suggest a combination of both. Com- 
paring Si+Au collisions at 14.5 A GeV and S+W collisions at 200 A GeV, both similar 
in size, we see a decrease in (/3f ) combined with an increase in the thermal freeze-out 
temperature (reflected in the flatter mj^-slope). This is connected with the increas- 
ing nuclear transparency at higher energies; the resulting fast longitudinal expansion 
leads to earlier freeze-out when transverse motion has not yet developed as strongly. 

Finally we consider the entropy per baryon S/B in the last row of Figure 16. 
There are two marks in each plot. The crosses correspond to the "measured" entropy 
obtained from the pion multiplicity^ or Dq. The triangles correspond to the values 
in a thermal model with the parameters plotted in the first 4 rows. The experimental 
values for S/B increase continuously from left to right in Figure 16. The strong rise of 
entropy production from AGS to SPS energies is expected. The observed size scaling 
can be understood in terms of increasing stopping power. We already discussed in 
Section 7, but want to point out again that the thermal model predictions lie above 
the measured values in Si+Au (AGS) collisions and N+N collisions, but below the 
measured values in the sulphur induced reactions at CERN. 

In summary, the thermal quantities can roughly be divided into two classes: Tf, 
Aq, As, and S/B scale dominantly with energy, while 7s and (/3f) scale mostly with 
size. We have seen that the qualitative behavior of these parameters is fully consistent 
with the expectations based on our present understanding of the heavy-ion collision 
dynamics. However, quantitatively the thermal hadron gas picture is plagued by 
internal consistencies. In the following Section we will pick up on these problems and 
use them to extract some further constraints on the early collision dynamics and the 
freeze-out process. 

9 Consequences for Freeze-Out Models 

In the preceding Sections 6-8 we presented an extensive study of (strange) parti- 
cle abundances in thermal models. Since a thermalized system surrounded only by 
vacuum necessarily begins after a while to expand hydrodynamically, and any ther- 
mal model analysis should for consistency be done on the basis of a hydrodynamical 
evolution background. This introduces the flow velocity as an additional essential pa- 
rameter into the analysis which must be determined consistently with the rest of the 
thermal parameters and with the dynamical evolution and freeze-out kinetics. In this 
light the strangeness signal for QGP formation acquires additional facets beyond the 
mere question of enhanced strange hadron production in ultrarelativistic heavy-ion 
collisions as originally discussed^^. As we showed, also the temperature and chemical 
potentials at the chemical freeze-out point can be extracted from strange particle 
abundances once the effects from the decay of unstable resonances after freeze-out 
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are properly taken into account. Combining the experimental information from the 
various strange particle abundances and from their momentum spectra we obtain 
valuable insight into the state of the system at freeze-out; using dynamical models 
we can then draw further conclusions also about earlier stages of the collision. 

We saw that, contrary to the situation at the AGS, at CERN energies the assump- 
tion of a thermalized hadronic system in relative chemical equilibrium cannot describe 
all observables. There is a crucial inconsistency with the observed pion production 
which we interpreted in terms of a lack of entropy in the hadron gas model compared 
to the data. Furthermore, the chemical analysis points towards chemical freeze-out 
temperatures near 200 MeV which are problematic for any weakly interacting hadron 
gas picture and at variance with lattice QCD results on thermally equilibrated sys- 
tems of strongly interacting particles. Up to now no fully consistent solution to these 
problems has been provided. However, there are a number of interesting theoretical 
suggestions in this connection which we would like to shortly discuss in this Section. 

9.1 Sequential Freeze-out 

An implicit assumption in the thermal and chemical analyses described in this review 
is that all particles freeze out simultaneously when the critical freeze-out temperature 
T = Tf is reached. This is a technical requirement related to the problem that 
it hard to include the modifications of partial freeze-out of some selected particle 
species on the subsequent hydrodynamical evolution of the remaining system. But 
it is nevertheless possible that different particle species decouple from the system at 
different times or temperatures, depending on their average cross sections with the 
medium-^-^'*'-^-^^ (which at CERN energies consists mainly of pions). Such a sequential 
freeze-out can occur both for the chemical composition of the fireball and for the 
momentum distributions. 

In Ref. 18 the idea of sequential freeze-out was applied to the measured strange 
and non-strange particle abundances. Freeze-out criteria of the type (103) lead to 
different freeze-out temperatures Tf for particles with different cross sections {av)[T). 
Cleymans et al.-^^ suggested that generically strange particles interact more weakly 
than non-strange hadrons and thus freeze out earlier. Assuming a typical factor 
2 between these two types of cross sections and adiabatic spherical expansion, they 
argue that freeze-out of strange particles at T = 185 MeV is consistent with freeze-out 
of non-strange particles at T = 130 MeV. The conceptual difference between chemical 
and thermal freeze-out temperatures was not considered in this work, and thus no 
theoretical justification for the observed similar mj^-slopes of all hadrons could be 
given. This model is able to resolve the excess entropy problem in terms of a larger 
specific entropy S/ B of the hadron gas at the lower freeze-out temperature of the non- 
strange hadrons; at the higher freeze-out temperature of the strange hadrons it agrees 
with the values around S/ B ~ 20 quoted above. It is not clear how this additional 
final state entropy arises in the model; implicitly it is generated between the strange 
particle and non-strange particle decoupling points. The question of viability of the 
hadron gas picture at the high strangeness freeze-out temperature is not addressed. 
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It is natural to ask whether there is any direct evidence for sequential freeze-out 
in nuclear collisions. It was suggested many years ago-^-^'*'-^-^^ that kaons should freeze 
out from a baryon rich environment earlier than most other hadrons, based on the 
very small K'^N cross section. As far as thermal freeze-out is concerned, this idea 
can be tested by looking at the two-particle momentum distributions ("Hanbury- 
Brown/Twiss (HBT) interferometry"): naively, if the kaons decouple earlier than, 
say, the pions and the fireball is expanding, the K^K^ correlation function should 
reflect a smaller HBT radius than the tt'^tt'^ correlation function. Measurements of 
kaon and pion correlations at the AGS^^^'^^^'^^^ and at the SPS^'^^'^'^'^ at first sight 
indeed seem to show such an effect. However, it is not clear whether this is indeed 
evidence for a smaller kaon freeze-out radius: Pion HBT radii are generically larger 
than kaon radii due to a larger contribution from delayed resonance decays which 
create a halo to the source. Furthermore, transverse collective expansion generates a 
characteristic l/y^mX dependence of the HBT radii^^^'^^^ which causes kaon radii to 
be smaller than pion radii if measured at the same p±. The situation at the AGS has 
not yet been fully analysed; at the SPS no statistically significant differences have 
been observed between K^K^ and K~ K~ correlations^^^ (for K~ the earlier freeze- 
out argument does not apply because it has a large cross section with nucleons), and 
the smaller kaon radius compared with the pion one is claimed to be consistent^^^ 
with the expected effects from resonance decays and transverse flow. 

Thus there is presently no clear direct experimental evidence for sequential thermal 
freeze-out. To test this idea further for chemical freeze-out requires extensive kinetic 
and hydrodynamic simulations. 

9.2 Sudden Hadronization Scenarios 

Another possibility to solve the inconsistencies of the thermal analysis is to postulate 
that that the particles do not decouple from an equilibrated system. However, even 
in that case one still has to explain the approximately thermal nature of the momen- 
tum spectra and the remaining order in the particle ratios which is described by the 
vanishing strange quark chemical potential, /ig = 0. A completely non-equilibrium 
approach to hadronization was recently initiated by the Budapest group^^^'^^°. Un- 
fortunately, this model cannot predict the final hadronic momentum spectra and is 
not yet in a form where it can be implemented into dynamical calculations. 

In Ref. 25 it was argued that the apparently large specific entropy and vanishing 
fis could be due to the existence of a QGP in the early stages of the collision which 
hadronizes suddenly into a hadronic system which immediately decouples. The exis- 
tence of a QGP would naturally explain the excess specific entropy and the van- 
ishing fis. Sudden hadronization followed by immediate decoupling was required in 
order to preserve the values of the chemical potentials across the phase transition and 
avoid their re-equilibration towards equilibrium hadron gas values before freeze-out. 
However, no consistent dynamical model which would produce the desired results was 
found in Ref. 25. 

The idea of sudden QGP hadronization (also called "QGP break-up model") was 
recently picked up by several groups^^^"^^^. They revived an old idea by van Hove^^'* 
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who suggested that, if the hadronization phase transition is of first order, hadroniza- 
tion could occur dynamically through a deflagration or detonation^ shock front. Such 
a violent process results in very rapid hadronization in which freely moving particles 
are ejected into the surrounding vacuum with large average velocities. 

The hadronization of a QGP through a shock front must satisfy energy-momentum 
and baryon number conservation across the front. This defines the Taub adiabat^^^ 

"•QGP V _ (eqcp + ^'QGp)(eQGP + ^h) ^^^^^ 



"•H / (en + PH)(eH + Pqgp) 

where n is the baryon number density, e the energy density and P the pressure in 
the initial QGP and final hadronic state, respectively. Even if for the final state an 
equilibrium hadron gas equation of state is assumed, Eq. (104) is not sufficient to 
determine the temperature and chemical potential change through the hadronization 
front. As a further constraint one often assumes entropy conservation which leads to 
the additional equation (Poisson adiabat) 

|^ = |i. (105) 

For shock fronts Eq. (105) is generally not a good assumption, but it can be used to 
give upper or lower limits on the changes of thermodynamic parameters during the 
hadronization. 

The idea behind the approaches in Refs. 121-123 is to see whether a detonation 
or deflagration scenario can provide a dynamical setting in which a superheated equi- 
librated hadron gas can arise in the final state which has T ~ 200 MeV, fis = 0, and 
is strangeness neutral, but which successfully fights its inherent tendency to return 
into the QGP state and rather decouples immediately due to its rapid dynamical 
expansion. However, these papers do not address the problem that such an equili- 
brated hadron gas would again yield too little specific entropy. In this sense they are 
complementary to the approaches of the previous subsection. No attempts have been 
made yet to include suitable deviations from thermodynamic equilibrium in the final 
state. 

Eqs. (104) and (105) define a detonation or deflagration front in space-time which 
can be either timelike-^^^ or spacelike-^^-^'-^^^. The time evolution is roughly as follows. 
The system starts from a QGP with a high initial temperature (assumed to be around 
200 MeV or higher). Then it supercools to about 0.6 — 0.8 Tc. Rapid longitudinal 
expansion and cooling of the QGP, combined with a very large nucleation time^^^ of 
around 100 fm/c for hadron gas bubbles, prevent the system from hadronization at 
Tc. The strong supercooling is basically caused by the requirement (105): only if the 
QGP is strongly supercooled and the hadron gas is strongly superheated the larger 
specific entropy in the QGP phase can be absorbed by the outgoing hadrons. 



shock front is called a detonation if in its rest frame the velocity of the incoming (QGP) 
matter is lower than that of the outgoing one (hadron gas), otherwise it is called a deflagration. 
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The supercooled QGP is mechanically unstable. Within the bag model, the QGP 
pressure turns negative^^^ already at T < 0.98 Tc. At 0.6 — 0.8 Tc sudden hadronization 
via a deflagration/detonation shock sets in. The outcome is a superheated hadron gas 
with a temperature around 200 MeV. After the deflagration/detonation the particles 
can considered as frozen out^^^'^^^. 

Models of this type explain naturally some of the puzzling experimental results: 

• fiB and fis are nearly conserved during the detonating hadronization^^^ . Thus 
the naturally vanishing strange quark chemical potential of a strangeness neutral 
QGP can be dynamically transferred to the final hadronic state. 

• Lattice QCD calculations'*^ favor a Tc around 150 MeV at /xb = 0. However, 
strangeness neutrality in the thermal hadron gas model with /ig = requires 
temperatures around Tf = 190 MeV. The mj^-slopes of all particles measured at 
CERN energies result in effective temperatures above 200 MeV. The superheat- 
ing to about 1.3 Tc in the calculation of Ref. 122 combined with the collective 
ejection of the hadrons from the shock front could explain all these features. 

• A QGP is a natural source of a large specific entropy^°°'^^°. Instead of (105), 
a realistic detonation or deflagration scenario would lead to further entropy 
increase^^'*. However, as already mentioned, this additional entropy can only 
be used profitably in comparison with the data if the final state is not modeled 
by an equilibrium hadron gas. In this respect all suggested models fail so far. 

• Assuming boost-invariant longitudinal expansion, HBT interferometry results^^^ 
suggest for sulphur induced collisions at 200 A GeV a total collision lifetime of 
about 4.5 - 6.5 fm/c. This agrees with the lifetime of the QGP phase until 
detonation from the supercooled state in Ref. 121. Assuming near-equilibrium 
hadronization by bubble nucleation would take much longer^^^. 

• The detonation into freely moving particles would naturally explain the high 
antiparticle yields seen in the experiments (see references in Table 9). Dur- 
ing slow hadronization the produced antiparticles would be absorbed in the 
longlived hadron gas subphase during and after the phase transition^ ^ . 

These points seem make the detonation/deflagration scenario an attractive can- 
didate for the dynamics of heavy-ion collisions at CERN energies. However, as far as 
we know, it only works if the hadronization is a strong first order phase transition. 
This is not supported by modern lattice QCD calculations, even when including finite 
net baryon number according to the best present knowledge'*^. Clearly much work 
is needed to supplement the present qualitative considerations by a fully consistent 
dynamical hadronization calculation. 
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10 Conclusions and Outlook 



We have reviewed new developments concerning the strangeness signal in heavy-ion 
experiments. In spite of tremendous experimental progress a clear proof of QGP for- 
mation still does not exist. However, circumstantial evidence in favor of a deconfined 
quark-gluon state in the initial stages of the collision, in particular at CERN ener- 
gies, keeps accumulating. A high degree of strangeness equilibration (75 ~ 0.7 — 1) 
indicates the opening of new, unusually rapid strangeness production channels. The 
mechanism could be QGP formation, but strong medium effects in a dense hadronic 
environment, leading to drastically reduced strange hadron masses and modified in- 
elastic cross sections, have also been shown to reproduce some of the observations. 
More detailed and comprehensive studies should be performed for further differenti- 
ation. 

The rather stable result of a vanishing strange quark chemical potential, fis = 0, 
in all of the CERN experiments further points in the direction of a state in which the 
phase space for strange quarks and antiquarks is symmetric. In any type of confined 
state this symmetry is spoiled by a finite net baryon density because the excess of 
light quarks over light antiquarks rubs off on the strange quarks because they are 
clustered together with them into hadrons. 

The observation of unexpectedly high pion production in heavy-ion collisions at 
CERN^^° points to a new efficient mechanism for entropy production as one passes 
from the lower BEVALAC and AGS energies to the SPS. 

We showed that all of the AGS data and most of the SPS data can be efficiently 
described by a thermal hadron gas model. The model contains only a few parameters 
with obvious intuitive interpretation. As shown in Section 8, they provide an efficient 
parametrization of the multitude of data, resulting in a comprehensive understanding 
of the basic features of the collision dynamics in ultrarelativistic heavy-ion collisions. 
The model can be tested with high accuracy using the large body of available data on 
particle yields and momentum spectra. This test is successfully passed by the data 
from Si+Al/Au collisions at the AGS, providing convincing evidence for the formation 
of a thermally and chemically equilibrated state with a large degree of strangeness 
saturation and showing clear signs of collective behavior. But even where the model 
fails in its application to hadron production data at CERN, it still provides a useful 
parametrization of the failure, thus leading to useful insights and stimulating further 
model building, as discussed in Section 9. 

The future will provide us with many new sets of data, especially from the larger 
collision systems Au+Au at the AGS and Pb+Pb at CERN. These will allow for fur- 
ther tests and confirmation of the picture developed so far. We expect even stronger 
evidence for equilibration and collective behavior at the AGS and perhaps a clarifi- 
cation of the remaining puzzles about the nature of the fireball at CERN. 

On the theoretical side, a major effort in constructing realistic dynamical mod- 
els for the hadronization of a QGP is required, in order to reach a comprehensive 
qualitative and quantitative understanding of the particle production process and to 
supplement the microscopic kinetic models presently under development^^ by a suc- 
cessful semiphenomenological macroscopic approach. Also the kinetics of the hadronic 
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freeze-out process and its proper implementation into macroscopic models should be 
studied in more detail. 

Finally, more efforts are needed to better understand the perturbative and non- 
perturbative contributions to strangeness production in a deconfined (equilibrium or 
non-equilibrium) quark-gluon state. 

Acknowledgments 

We would like to thank Marek Gazdzicki, Jean Letessier, Ahmed Tounsi and 
especially Johann Rafelski for several years of collaboration on this subject. Many of 
new results reported in this overview were obtained by them or in collaboration with 
them. This work was supported by DFG, BMBF, and GSI. 

A Appendix: Literature Guide to Strange Particle Mea- 
surements 

Although this article is meant as a theoretical review on strangeness, the phe- 
nomenology of strangeness requires the input from and the comparison with experi- 
ment. For this reason we present here a list of experimental publications on strange 
particle production in nuclear collisions, organized by particle species and experiment 
number. The two tables summarize the experiments at the Brookhaven AGS and the 
CERN SPS, respectively. We concentrate on strange hadron spectra and multiplicities 
and left out articles on e.g. kaon interferometry. 

The list mainly contains articles in refereed journals. Conference proceedings are 
mostly neglected with the exception of Quark Matter '91"^^^ and '93"^^^. We recom- 
mend for updated strangeness data also the proceedings of Quark Matter '95^^ and 
Strangeness '95^^, of which a few contributions have already been included according 
the notes taken by one of the authors (UH) and which will be published soon. 

The tables should contain all relevant publications from the last few years known 
to the authors; it is, however, quite likely that we missed one or the other important 
contribution for which we would like to apologize. 
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